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Abstract 

A fundamental problem in computational geometry is to compute an obstacle-avoiding Eu- 
clidean shortest path between two points in the plane. The case of this problem on polygonal 
obstacles is well studied. In this paper, we consider the problem version on curved obstacles, 
which are commonly modeled as splinegons. A splinegon can be viewed as replacing each edge 
. of a polygon by a convex curved edge (polygons are special splinegons). Each curved edge is 

04 ' assumed to be of 0(1) complexity. Given in the plane two points s and t and h splinegons of 

5-j . totally n vertices with pairwise disjoint interior, we compute a shortest s-io-t path avoiding the 

splinegons, in 0{n -\- hlo^'^'^ h + k) time for any arbitrarily small constant e > 0, where fc is a 
parameter sensitive to the structures of the input splinegons and is upper-bounded by 0{h^). In 
particular, when all splinegons are convex, k is proportional to the number of common tangents 
^Sj . in the free space (called "free common tangents") among the splinegons. Previous work is known 

only for some special cases of this problem (e.g., when the obstacles are discs or certain convex 
objects). We develop techniques for solving the problem on the general splinegon domain which 
also improve several previous results. (1) We improve the previously best solution for computing 
\^ . a shortest path between two points among convex pseudodisks of 0(1) complexity each, which 

t/3 ' find applications in robot motion planning. (2) We even improve an algorithm for the polygon 

case (i.e., when all splinegons are polygons). The polygon case has been solved in O(nlogn) 
time (the worst case optimal), or in 0{n + h^ logn) time (this was also for the special case when 
all polygonal obstacles are convex). Thus, we improve the 0{n + logn) time result. (3) We 
give an optimal algorithm for a basic visibility problem of computing all free common tangents 
among h pairwise disjoint convex splinegons of totally n vertices, in 0{n + k + hlogh) time 
' and 0{n) space, where k is the number of all free common tangents (and k = 0{h'^)). Even for 

fT^ , the special case when all splinegons are convex polygons^ the previously best visibility algorithm 

' takes 0(n -I- /i^logn) time. We also present an optimal algorithm for computing the Voronoi 

, diagram of pairwise disjoint convex splinegons. Our techniques may find other applications. 

1 Introduction 
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Finding an obstacle avoiding Euclidean shortest path in the plane between two points is a funda- 
mental problem in computational geometry and has many applications. The case of this problem 
on polygonal obstacles is well studied (e.g., [HI [M 13 Ell [Ml EH IS] ) • For obstacles bounded by 
curves, the problem is more difficult and only limited work is found in the literature. In this paper, 
we study the problem version on curved obstacles and present efficient algorithms for it. 

1.1 The Geometric Setting and Our Results 

As in [71 |8l [27], we use splinegons to model planar curved objects. A (simple) splinegon 5 is a 
simple region formed by replacing each edge e- of a simple polygon P by a curved edge Cj joining 
the endpoints of e- such that the area bounded by the curve Cj and the line segment is convex 
(see Fig. [T|). The vertices of S are the vertices of P. Normally, each curved edge of S is represented 
parametrically. Note that computation on curved objects is a challenging research topic (e.g., see 
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Figure 1: A splinegon (solid curves) defined on a poly- Figure 2: The relevant visibility graph of three convex 
gon (red or dashed segments). objects. 



[3]). In this paper, as in [51 127], we assume that each sphnegon edge is of 0(1) complexity, and 
primitive operations on a sphnegon edge can each be performed theoretically in 0(1) time (e.g., 
computing its intersection with a line). 

We study the problem of computing shortest paths in a splinegon domain, denoted by SPSD. 
The plane minus the interior of obstacles is called the free space. Given two points s and t and a 
set of h splinegons, S = {Si, . . . , Sh}, of totally n vertices with pairwise disjoint interior, the SPSD 
problem seeks a shortest path from s to t in the free space. If the splinegons in S are all convex, 
then we refer to this case as the convex SPSD. To our best knowledge, we are not aware of any 
previous work for the general SPSD. On the convex SPSD, by generalizing the algorithm in |32) 
for the convex polygonal domain, one may be able to obtain an 0{n + h"^ logn) time solution. 

We develop techniques for the general SPSD problem and our algorithm, denoted by ^4/^0- 
SPSD, takes 0{n + /ilog^^"^ h + k) time, where A: is a parameter sensitive to the structures of the 
input splinegons and is bounded by 0{h'^) (the exact definition of k will be given in Section [2]). 
Throughout this paper, we let e > denote an arbitrarily small constant. For the convex SPSD, 
Algo-SPSD takes 0(n-|-/ilog/i-|-A:) time, with k being the number of free common tangents among 
the splinegons. A common tangent of two convex splinegons is a line segment that is tangent to 
both splinegons at its endpoints; the common tangent is free if it lies entirely in the free space. 

We also solve a relevant visibility graph problem: When all splinegons in S are convex, compute 
the free common tangents of the splinegons (see Fig. [2]). Our algorithm for it takes optimally 
0{n+h\og h+k) time and 0(n) space. This visibility problem is a key subproblem to our algorithm 
Algo-SPSD, and is interesting in its own right. In addition, we present an optimal algorithm for 
computing the Voronoi diagram of a set of pairwise disjoint convex splinegons. 

1.2 Related Work 

The polygon case of SPSD (i.e., S contains polygons only) is well studied. By constructing the 
visibility graph |12] . a shortest s-t path can be found in 0(nlogn + K) time, where K = 0{v?) 
is the size of the visibility graph. By building a shortest path map, Storer and Reif solved this 
case in 0{nh) time [33]. Mitchell |28] gave the first subquadratic, 0(n^/^"'"'^) time algorithm for 
it based on the continuous Dijkstra approach. Also using the continuous Dijkstra approach and 
a conforming planar subdivision, Hershberger and Suri |16| presented an O(nlogn) time solution. 
An 0{n + /i^ log n) time algorithm was given in [21] (a preliminary version is in [19] and full details 
are in [20]). Thus, our Algo-SPSD algorithm improves the results in |2H 133] and is faster than the 
O(nlogn) time solution [16] for small value h, say h = o{y/nlogn).* 

For SPSD on curved obstacles, only limited results for some special cases are known. For the 
case with n discs, O(n^logn) time algorithms were given [31 16], and a heuristic approach [23] was 
derived with experimental results. For disks of the same radius, the algorithms in ^5] 133] can find a 
shortest s-t path in 0{in?) time. A set of objects in the plane is called pseudodisks if the boundaries 

*Very recently, an unrefereed report ^8J announced an algorithm for the polygon case based on the continuous 
Dijkstra approach with an 0{n + /i log ft log n) time. Our algorithm is superior to it when k = o(n + ft log /i log n) . 



2 



of any two objects can cross each other at most twice. If S contains n convex pseudodisks of 0(1) 
complexity each, an algorithm in [5] can find a shortest s-t path in 0{n?) time. By using om' 
Algo-SPSD algorithm, the result in [5] can be improved as follows. Let iyJ5 denote the union of 
the convex pseudodisks in S and K be the number of vertices on the boundary of y5. It has been 
shown in [22] that K = 0{n) and iyJ5 can be computed in 0(n log^ n) time. Since all pseudodisks in 
S are convex, iyJ5 can be viewed as consisting of pairwise disjoint splinegons; let H be the number 
of splinegons in iyJ5 (obviously, H < n). By applying Algo-SPSD to iyj5, a shortest s-t path can 
be found in 0(n log^n + k) time, where k is bounded by 0{H^). This improves the O(n^) time 
result in [5] when H = o{n). As a consequence, a robot motion planning problem [Sj [15] can also 
be solved in a faster manner. 

For a single splinegon S, a shortest s-t path in S can be found in 0{n) time, and further, 
shortest paths from s to all vertices of S can be found in 0{n) time [27| . 

It is not clear how to apply the continuous Dijkstra approach [161 128] to our SPSD problem 
(even when all splinegons are discs) due to the curved obstacle boundaries. For example, Mitchell's 
approach [2^ uses a data structure for processing wavelet dragging queries by modeling them as 
high-dimensional radical-free semialgebraic range queries. In SPSD, however, such queries would 
involve not only radical numbers but also inverse trigonometric operations (e.g., arcsine), and 
hence similar techniques do not seem to apply. Hershberger and Suri's approach |16j relies heavily 
on a conforming subdivision defined by the vertices of the polygonal obstacles. In our problem, 
however, it seems highly elusive to determine a set of 0{n) vertices or points that can help build 
such a subdivision. One might attempt to use the splinegon vertices to build such a subdivision. 
But, an important property used by the subdivision [16] is that the generator of every wavelet must 
be one of the obstacle vertices. Yet in SPSD, a generator need not be a splinegon vertex. 

Constructing the visibility graph for polygonal objects is well studied [T| W2 \ [TH l [1^ 1 [52 } [5^. 
The fastest algorithm for it takes 0(n log n -|- K) time [12], where K = 0{'n?) is the size of the 
visibility graph. For the relevant visibility graph problem [3H [32] (or building the relevant 
visibility graph) on splinegons, two special cases have been studied. When S contains n convex 
objects of 0(1) complexity each, the problem is solvable in 0(nlogn -|- K) time [31], where K = 
0{n'^) is the number of free common tangents. If S contains h convex polygons, as in [2H I32|. then 
the problem is solvable in 0{n + /i^logn) time; an open question was posed in [21j to solve this 
case in 0{n + k\ogn) time, where k = 0{h'^) is the number of free common tangents. Note that 
our optimal 0{n + h log h + k) time result is better than the solution desired by this open question. 

2 An Overview of Our SPSD Algorithm 

Our algorithm Algo-SPSD follows the high level scheme used in the polygonal domain case [21] , 
but with the key steps replaced by our new, generalized, and more efficient solutions for the more 
difficult splinegon domain counterparts. Since the paper is long, we present it in a way that each 
section is highly self-contained, as discussed below. Let 7^ be a rectangle containing all splinegons 
in 5, and J- denote the free space inside TZ. We view both s and t as two special splinegons in S. 

The first step is to decompose J- into regions each of which has at most four sides and at most 
three neighbors (see Fig. [3]). This decomposition, called bounded degree decomposition, serves the 
same purpose as a usual triangulation in the polygonal domain case. Melissaratos and Souvaine [27] 
computed a bounded degree decomposition inside a simple splinegon in linear time. By extending 
the triangulation algorithm for the polygonal domain case [2j and applying the algorithm in [27J, 
we present an 0{n + /ilog^"*"*^ h) time algorithm for computing a bounded degree decomposition of 
F, denoted by BDD[F). The details of this step are given in Appendix [Cl 
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Figure 3: Illustrating a bounded degree decomposition of T Figure 4: Illustrating an open hourglass (left) and 

(with dashed segments) and the corridors (with red solid arcs), a closed hourglass (right) with a corridor path linking 

There are two junction regions indicated by large (red) points the apices x and y of the two funnels. The dashed 

inside them, connected by three solid (red) arcs. Removal of segments are diagonals. The paths 7r(a, fe) and7r(c, d) 

these two junction regions results in three corridors. are shown with thick solid curves. 

The second step, with details in Appendix [Dl is to compute a corridor structure in BDD[T). 
The corridor structure consists of 0{h) corridors and 0{h) junction regions (see Fig. [3]). Each 
corridor contains an hourglass, either open or closed (see Fig. [4|). An open hourglass contains two 
convex chains. A closed hourglass contains two "funnels" with a corridor path connecting the two 
apices of the two funnels. Each side of a funnel is also a convex chain. As in [2T], the above 0{h) 
convex chains from the corridors can be used to partition into a set S' of 0{h) convex splinegons 
of totally 0{n) vertices with pairwise disjoint interior such that a shortest s-t path for our original 
SPSD problem is also a shortest s-t path avoiding the interior of the convex splinegons of S' and 
possibly utilizing some corridor paths. Thus, in addition to the presence of the 0{h) corridor paths, 
our SPSD problem is reduced to an instance of the convex SPSD. All the above computation can 
be performed in 0{n + /ilog h) time. The key is to solve the convex SPSD problem on S' . 

To solve the convex SPSD on S' , we define a relevant visibility graph G (see Fig. [2|), as follows. 
Let k be the number of all free common tangents of the 0{h) convex splinegons in S'\ thus k = 
0{h'^). The node set of G consists of the endpoints of the free common tangents. Hence G has 
0{k) nodes. Each free common tangent defines an edge in G. For every splinegon S € S\ its 
boundary portion between any two consecutive nodes of G along the boundary of S also defines an 
edge. Thus G has 0{k) edges. Clearly, a shortest s-t path in the free space of S' corresponds to a 
shortest path from s to t in G (both s and t are nodes in G). Therefore, to solve the convex SPSD, 
we need to solve two subproblems: Constructing G and computing a shortest s-io-t path in G. 

The third step solves the first subproblem: Building G. Our algorithm takes 0{n + k + h log h) 
time. This step, the most challenging part of our SPSD algorithm, is presented in Section [3j 

The fourth step solves the second subproblem: Finding a shortest path from s to t in G. Since 
G has 0{k) nodes and 0{k) edges, simply running Dijkstra's algorithm on G would take 0(A;logA;) 
time. To avoid the log k factor, we extend the approach in [5j for computing a shortest path among 
pseudodisks. A problem we need to handle is to compute the Voronoi diagram of the convex 
splinegons in S' . We show that this Voronoi diagram can be computed in 0{n + hlogh) time, 
which may be of independent interest. The details of this step are discussed in Appendix lEl 

The last step, discussed in Appendix |Fl is to incorporate the 0(h) corridor paths into the 
algorithm for the convex SPSD on S' to obtain a shortest path for our original SPSD problem. 

3 Computing the Relevant Visibility Graph of Convex Splinegons 

In this section, we construct the relevant visibility graph G for the 0{h) convex splinegons in S' of 
totally 0{n) vertices. For convenience, we slightly change the notation and consider the following 
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Figure 5: A flip operation on a free bi- Figure 6: The common Figure 7: R = Rtri{b) and h = ha. 

tangent 6: (a) The dashed bitangent is b tangent hne of two pseudo- The awake points are those on x^qii. 

before the flip; (b) the dashed bitangent triangles. To flip 6, the walk on dR starts at q'n 

is if{b) after the flip. for Case 1 and at xr for Cases 2 and 3. 

problem: Construct the relevant visibility graph G for h pairwise disjoint convex splinegons in a 
set V = {Pi, P2, ■ ■ ■ , Ph} of totally n vertices. Let B denote the set of all free common tangents of 
V and let k = \B\. Our algorithm for computing B runs in 0{n + k + h log h) time and 0{n) space. 
The graph G can also be constructed in 0(n + k + h\og h) time. 

Our algorithm can be viewed as a generalization of Pocchiola and Vegter's algorithm [31j (and its 
preliminary version [30j). We call it the PV algorithm. Given n pairwise disjoint convex obstacles of 
0(1) complexity each, O = {Oi, O2, ■ ■ ■ , On}, the PV algorithm computes all free common tangents 
of O in 0(nlogn + K) time and 0(n) space, where K = 0{in?) is the number of all free common 
tangents of the n obstacles in O. It was also claimed in [31] (without giving any details) that the 
PV algorithm may be made to compute B for our problem on V in 0(n log h + k) time. 

In general, the PV algorithm relies mainly on the convexity of the obstacles involved. The 
needed properties also hold for our setting. The high level scheme of our algorithm follows that 
of the PV algorithm, but with certain modifications. However, the most challenging task is to 
achieve a tight and optimal time bound for the algorithm. For the problem in [3T| solved by the 
PV algorithm, this is easy. For our problem, the same analysis as for the PV algorithm does not 
work. As shown later, our analysis needs numerous non-trivial new observations and novel ideas. 
Comparing with the PV algorithm, our algorithm and analysis explore more crucial properties and 
structures of the problem, which may be useful for solving other related problems as well. Further, 
our modifications to the scheme of the PV algorithm seem necessary (without them, it is not clear 
to us whether the scheme in [31] can be generalized to attain our optimal time bound). 

3.1 An Overview of Our Relevant Visibility Graph Algorithm 

We maintain a structure, called "good pseudo-triangulation" , as a subdivision of the free space by 
some noncrossing free common tangents (e.g., Fig.[5|). After an initial good pseudo-triangulation is 
computed, we repeatedly choose a certain free common tangent b in the good pseudo-triangulation 
and "flip" b (somewhat related to a flip on a triangulation) to obtain a new free common tangent 
b* and a new good pseudo-triangulation. After k flips, all common tangents in B are computed. 
The order of the flips is guided by a topological (partial order) structure maintained in these good 
pseudo-triangulations. The effectiveness of this approach hinges on its ability to perform the k flips 
in 0{n + k) time, and the key is to determine b* efficiently. The main idea is discussed below. 

A good pseudo-triangulation contains many pseudo-triangles. For each pseudo-triangle T, we 
choose and store a critical portion of its boundary, as Awake\T], which is used to find b* . After 
obtaining b* and a new good pseudo-triangulation, we need to update Awake for some (two) new 
pseudo-triangles induced by the flip. To update Awake efficiently, we also choose and store a 
boundary portion of each pseudo-triangle T, as Asleep[T]. In other words, Awake is used to find 
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b* and Asleep is used to update Awake {Asleep itself also needs to be updated). A key difference 
between the PV algorithm and ours is that Asleep[T] refers to different portions of T's boundary. 
Without this modification, it is not clear to us how to generalize the original PV algorithm. 

Both Awake and Asleep are implemented as "splittable queues" [H] that support three opera- 
tions: Enqueue, dequeue, and split. Our algorithm uses two phases for handling each flip: Phase I 
computes b*; Phase II updates Awake and Asleep. To bound the overall running time, it suffices to 
prove the following key claim: The total number of enqueue operations for all k flips is 0(n + k). 
Actually, for each flip, only 0(1) sequences of enqueue operations are needed and each sequence is 
on either a free common tangent or a boundary portion of a single obstacle. 

A most essential difference between the PV algorithm (on the problem in [31]) and ours is 
on proving the key claim. In the PV algorithm, this is fairly easy: Since every obstacle is of 
0(1) complexity, each enqueue sequence needs only 0(1) enqueue operations. In our problem, 
however, this is quite challenging. This is because the complexity of the boundary of an obstacle 
(i.e., a splinegon) can be Q(n), and thus an enqueue sequence may take as many as 0(n) enqueue 
operations. To prove the key claim, we must analyze it in a global manner, which needs many new 
observations and analysis techniques and is one of our main contributions in this paper. Appendix 
IB] is devoted entirely to this task. 

3.2 Terminology and Basic Geometric Structures 

We focus on showing how to compute B since the graph G can be constructed accordingly while 
B is being computed. In the following, we simply call each splinegon in V an obstacle and each 
(curved) splinegon edge an elementary curve. Thus, the complexity of each elementary curve is 
0(1). We call a common tangent of two obstacles a bitangent. As in the PV algorithm [30\ I31j . 
we assume that all obstacles in V are smooth (i.e., only one tangent line touches each boundary 
point) and are in general position (i.e., no three obstacles share a common tangent line). These 
assumptions are only for ease of exposition. To handle the case with polygons, for example, we can 
take the Minkowski sum of the polygons with an infinitesimally small circle. 

As in [30^ I31|. we define a pseudo-triangulation of the convex obstacles in "P as a subdivision of 
the free space induced by a maximal number of pairwise noncrossing free bitangents (see Fig. [5]). 
As shown in [30\ I31j . the number of free bitangents in any pseudo-triangulation of V is 3h — 3. 

Let 7~ be a pseudo-triangulation and B(T) denote the set of all free bitangents that appear in 
T. Any bounded free face T in T is a pseudo-triangle, and the boundary of T, denoted by dT, 
consists of three convex chains with convexity towards the interior of T. The three endpoints of 
the convex chains are called the cusps of T. Denote by B{dT) the set of all free bitangents on dT 
(i.e., each bitangent in B{dT) lies on dT). For a point p on dT lying on an obstacle Pj, the tangent 
line of Pi at p is called the tangent line of T at p; for a point p lying on a free bitangent of B{dT), 
the line containing the bitangent is the tangent line of T at p. Note that if a line is tangent to 
T at a point p G dT, then the line is also tangent to the convex chain of dT that contains p. As 
shown in \30\ I31j . any two pseudo-triangles in T have a unique common tangent line, i.e., a line 
tangent to both pseudo-triangles (see Fig. [6|). Suppose two adjacent pseudo-triangles T and T' in 
T share a bitangent b G B{T); a flip operation on b replaces b by the common tangent of T and 
T' , which is a free bitangent and denoted by 99(6) (see Fig. [5]). A flip operation produces another 
pseudo-triangulation. If b lies on the convex hull of V, then we let '.p{b) be b itself. 

Our algorithm for computing B, called the topological flip algorithm, performs flip operations 
based on a topological order, which can be viewed as a generalization of a topological sweep |9j. 
We define a partial order on the free bitangents and a topological structure (like the horizon trees 
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in [9]) which is maintained by our algorithm. To define the topological structure, below we assign 
directions to bitangents and tangent lines of pseudo-triangulations, and discuss some properties. 

Given a unit vector u, the u-slope of a directed line (or segment) / is defined as the angle 
(in [0,27r)) of rotating u counterclockwise to the same direction as /. For an undirected line (or 
segment) /, its u-slope is the angle (in [0,7r)) of rotating u counterclockwise to the first vector 
parallel to I; the direction of that vector is said to be consistent with the n-slope of /. 

Consider a pseudo-triangulation T. Given a vector n, for each bitangent b S B{T), we assign 
to b the direction consistent with the u-slope of b. For every pseudo-triangle T of T, let br be 
the bitangent in B[dT) with the minimum u-slope. Further, for each point p G 5T, we assign to 
the tangent line l{p) of T at p the direction consistent with the ^y-slope of l{p), and call l{p) the 
directed tangent line of T at p. The 6T-slope of the directed l{p) is also called the pseudo-triangle 
slope (or pt-slope for short) of l{p) at the point p £ dT. For any bitangent b € B, suppose we assign 
a direction to b and p is an endpoint of b (say, p is on dT of a pseudo-triangle T); then the direction 
assigned to b is said to be compatible with p if the directed tangent line of T at p has the same 
direction as b. Note that the pt-slope of any point on 6t is zero. As moving on dT clockwise from 
6t, the pt-slope of the moving point increases continuously from to it, until we are back to bT- 

Our algorithm maintains a topological structure called good pseudo-triangulation, defined as 
follows. A pseudo-triangulation T is said to be good (called weakly greedy in [30]) if there is a way 
to assign every free bitangent b G B a. direction such that a partial order -< can be defined on the 
directed bitangents of B{T) with the following properties: (1) For each pseudo-triangle T in T, 
the partial order ^ is a total order, which corresponds to the pt-slope order on B{dT) with respect 
to 6t; (2) the direction of each bitangent b G B \s compatible with both its endpoints; (3) for any 
bitangent b G B\ B{'T), all bitangents in B{T) intersecting b cross the directed b from left to right. 

3.3 Initialization and Processing of the Topological Flip Algorithm 

Let uq be a vector with the direction of the positive x-axis. We first compute an initial pseudo- 
triangulation To of P induced by a set {61,62, . . . ,63/1-3} of free (undirected) bitangents such that 
(1) 61 is the bitangent in B with the smallest uo-slope, and (2) for any 1 < i < 3/i — 3, 6j_|_i is the 
bitangent with the smallest no-slope in B that does not cross any of 61,62, ■ ■ ■ ,bi (e.g.. Fig. [5|^a) is 
To)- As shown in [301 [31], To for the obstacle set O is a good pseudo-triangulation, and can be built 
in 0(n log n) time. Likewise, for our problem. To of "P is also a good pseudo-triangulation; further, 
we can compute To even faster, as shown in Lemma [1] (the proof is in Appendix lUj) . 

Lemma 1 The initial good pseudo-triangulation TqoJV can be constructed in 0{n-\- hlogh) time. 

After computing To, we assign to every bitangent 6 G B(To) the direction consistent with its 
uo-slope (in [0,7r)). For each bitangent b £ B \ B{To), as shown in [3Ql|31], we can always assign a 
direction to 6 that is compatible with both 6's endpoints; for the purpose of discussion, we assume 
that this direction has been assigned to 6 (actually, the algorithm does not perform this assignment 
explicitly). Let T = To- To compute the bitangents in B\B{To), the topological flip algorithm keeps 
flipping a minimal bitangent in B{T) and generating another good pseudo-triangulation of 7^. A 
(directed) bitangent 6 G B{T) is minimal \i it has the smallest tio-slope among all free bitangents on 
the boundaries of both the left and right adjacent pseudo-triangles of 6 in T. A minimal bitangent 
always exists in a good pseudo-triangulation T [30l [Si]- The next lemma (proved in [30] and 
applicable to our problem) shows that any minimal bitangent in B{T) can be flipped. 

Lemma 2 jlSOf For any good pseudo-triangulation T of V (initially, T = To), let b be a minimal 
free bitangent in B{T) with a uq- slope less than vr. Then the pseudo-triangulation T' ofV obtained 
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by flipping b is also good, and the assigned direction of any free bitangent t € B \ {b} does not 
change after the flip of b, whereas the direction of b is reversed. 

As shown in [301131). we only need to flip a minimal bitangent in B{T) with a tiQ-slope less than 
TT, and this ensures that the algorithm will terminate and all free bitangents in B will be generated. 
Note that once a bitangent is flipped, since its direction is reversed, its uo-slope becomes no smaller 
than vr, and thus it will never be flipped again. 

3.4 Conducting the Flips 

Given To, we determine the set of minimal bitangents in B{Tq), denoted by C. Then, we take an 
arbitrary bitangent b from C, flip 6, and update C. We repeat this process until C = (by then B 
is obtained). The key is to perform all k (= \B\) flips in 0(n + k) time. 

Let T be a good pseudo-triangulation. For any bitangent t in B{T), denote by Ltri{t) (resp., 
Rtri{t)) the pseudo-triangle of T (if any) that is bounded by the directed t and is on the left 
(resp., right) of t. Suppose we are about to flip a minimal bitangent b in T. Let R = Rtri{b) and 
L = Ltri{b). To compute (/?(6), an easy way is to walk clockwise along dR and dL synchronously, 
starting from 5, until finding (/7(6). But, this is too expensive. A more efficient approach is to first 
"jump" to a certain location on dR and dL and then do the synchronous walking. To implement 
this idea, we need some "crucial points" on dT for each pseudo-triangle T G T, as defined below. 

For any directed free bitangent 5, we denote its two endpoints by Tail{b) and Head{b), respec- 
tively, such that 6's direction is from Tail{b) to Head{b), and call them tail and head of b. 

Consider a pseudo-triangle T £ T. We define the basepoint of T, denoted by pT, to be the tail 
of bx (i.e., the smallest uo-slope bitangent in B{dT)) if T = Rtri^bj-), and be the head of 6^ if 
T = Ltri{bT)- Starting at pT, if we move along dT clockwise, the successive cusps of T encountered 
are denoted by xt, Vt, and zt (if pr is a cusp, we let it be zt)- The forward (resp., backward) 
T-view of any point p on dT is the intersection point of dT with the directed tangent line l{p) of 
T at p, lying ahead of p (resp., behind p). Let qx denote the special point on dT whose forward 
(resp., backward) T-view is px ^^T = Rtri{bT) (resp., T = Ltri{bT)) (e.g., see Fig. [7]). For any two 
points pi and p2 on dT, let piP2 denote the portion of dT from pi clockwise to p2- 

For a pseudo-triangle T, a point p G dT is said to be awake if and only if p G xtOt (see Fig. [7]). 
We let Awake[T] represent the awake portion xtQt of dT. Also, we let Asleep[T] represent the 
portion wtPt of dT, where the point wt = Qt if Qt £ Vtzt and wt = Vt if Qt ^ Vtzt- Note that 
our definition of Asleep[T] is different from that in [301 131]. in which Asleep[T] represents wtzt- 

To perform a flip on b, we use Awake[T] to find ip{b) and use Asleep[T] to help update Awake[T]. 
After every fiip, Awake[T] and Asleep[T] are updated accordingly. Both Awake[T] and Asleep[T] 
are stored as splittable queues, which support three types of operations on a list: (1) Enqueue an 
atom, either at the head or the tail of the list; (2) dequeue the head or the tail of the list; (3) 
split the list at an atom x, which is preceded by a search for x in the list. A data structure for 
implementing a splittable queue was given in [301 131| . whose performance is shown below. 

Lemma 3 ]3C^ A sequence of 0{n + k) enqueue, dequeue, and split operations on a collection 
of n initially empty splittable queues can be performed in 0{n -\- k) time. 

An atom in our problem can be a bitangent or an elementary curve. Further, an elementary 
curve may be divided into multiple pieces by the endpoints of some free bitangents in a good 
pseudo-triangulation, in which case an atom may be only a portion of an elementary curve. Li any 
case, the complexity of each atom is 0(1). For a pseudo-triangle T, a portion of dT is called an 
obstacle arc if that portion lies entirely on the boundary of a certain obstacle. 
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Figure 8: The three possible cases when flipping b; b* = Lp{b) and b'n is the next bitangent of b along dR clockwise. 

Initially, we compute Awake[T] and Asleep[T] for each pseudo-triangle T in To, using 0(n) en- 
queue operations. Consider a flip on a minimal bitangent b in the current good pseudo-triangulation 
T. Let R = Rtri{h) and L = Ltri{b). Let b* = ip{b), p* = Tail{b*), and q* = Head{b*). Let V be 
the resulting pseudo-triangulation after the flip. Let R' = Rtri[b*) and L' = Ltri{b*). To maintain 
the minimal bitangents in T' , we need to find the bitangent with the smallest lig-slope in B{dT) 
(i.e., for each T G Below, we only discuss the case for R' (the case for L' is similar). 

Let 6^ be the next bitangent of b in B{dR) clockwise along dR (i.e., b'j^ is the bitangent in 
B{dR) \ {6} with the minimum UQ-slope). Let be the portion of dR from Head{b) clockwise to 
the first encountered point of b'^. Clearly, is an obstacle arc, and b^i is one of the two bitangents 
6^ and b* (the one with a smaller UQ-slope). As in [3^ I31j. there are three main cases (see Fig. [5]). 

Case 1: b and 6^ are not separated by the cusp xr of R (i.e., Fr does not contain xji). Then 
R' (= Rtri{b*)) is also Rtri{b'^ and p* does not lie on F/j. Thus, bn' is b'j^. 

Case 2: b and 6^ are separated by xr and p* does not lie on Tr. Then 5^/ is b'^^. In this case, 

is either Head{b'^ (see Fig. [8]) or Head{b). As in [30l[31] and shown later, R' is also Ltri{b'^. 

Case 3: b and 6^ are separated by xr and lies on F^^. Then, 6^/ is 5* (e.g., see Fig. [8|). 

In the following, we divide the processing of a flip on b into two phases as in [30^ [3T] : Phase I 
finds b*; Phase II updates Awake and Asleep accordingly (i.e., computing AwakelR'], Awake[L'], 
Asleep[R'], and Asleep[L']). We only discuss the case for R' (the case for L' is similar). 

3.5 Phase I of a Flip Operation 

Let q'j^ be the point on dR whose forward i?-view is Tail{b'pj) (see Fig. [Tj); similarly, let be the 
point on dL whose backward L-view is Head{b'j^) (similar to 6^, b'^ is the next bitangent of b in 
B{dL) clockwise along dL). To compute b* = ^{b), Cases 2 and 3 are handled in the same way 
while Case 1 is different. For now, we assume that we already know whether it is Case 1, and 
further, when Case 1 occurs, we know b'^. We will show later how to detect the cases and find b'j^ 
for Case 1. An easy but useful observation is that p* = Tail{b*) must be on xrqr (= Awake[R\), 
and in Case 1, p* is even on q'^QR- Similarly, q* = Head{b*) lies on xiql (= Awake[L\), and even 
on q'l^qi in Case 1. To compute 6*, we walk clockwise along dR and dL synchronously, as follows. 
In Cases 2 and 3, the walk on dR starts at xr; in Case 1, we first split Awake[R] at q'p^ and then 
start walking from g^. The split operation on Awake\R\ at q'^ in Case 1 is preceded by a search for 
q'^ in Awake[R\, which is guided by the position of Tail{b'p) with respect to the directed tangent 
lines of R at the endpoints of the atoms in Awake[R\. Similar things are done on dL. 

We perform the following three main steps for computing the endpoints p* and q* of b* . 

Step (1): For Case 1, we split Awake[R\ (= xriir) at q'j^ into AwakeMin[R\ and AwakeMax[R\ 
such that the atoms in the former queue have the smaller pt-slopes; for Cases 2 and 3, we simply 
let AwakeMin[R] = and AwakeMax[R] = Awake[R]. Thus, in any case, p* always lies on an 
atom in AwakeM ax[R\. We compute AwakeMax[L] and AwakeMin[L] similarly. 

Step (2): Compute b*. To do so, the synchronous walks on dR and dL can be implemented by 
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dequeuing atoms from AwakeMax[R] and AwakeMax[L], until p* and q* are found. 

Step (3): Note that the atom at one end of AwakeMax[R] now contains p*. We cut that atom 
at p* by setting p* as an endpoint such that AwakeMax[R] now represents the portion p*qR. We 
do similar things for AwakeMax[L] and q* {AwakeMax[L] thus represents the portion q*qL)- 

These steps can be done at the cost of at most two split operations (at Awake[R] and Awake[L]), 
followed by multiple successive dequeue operations, but no enqueue operation is needed at all. 

To determine whether Case 1 occurs and find b'j^ in Case 1, we use the following preparing 
procedure. We walk from Head{h) along dR clockwise until first encountering either b'^ or xr. If 
we encounter 6^ first, then it is Case 1 and we have 6^ as well. Otherwise, it is Case 2 or Case 3. 

This finishes the description of Phase I. The discussion of Phase II is given in Appendix [Al 

3.6 Proving the Complexities of the Topological Flip Algorithm 

For ease of discussion, suppose we already have the following key lemma. 

Lemma 4 The total number o/ enqueue operations in Phase II of the entire algorithm is 0{n + k). 

Recall that the initialization procedure performs 0{n) enqueue operations and no enqueue 
occurs in Phase I at all. By Lemma [H the total number of enqueues in the entire algorithm is 
0(n + k). Thus, the total number of dequeues in the entire algorithm is also 0{n + k) since it 
cannot be bigger than the total number of enqueues. Further, at most two splits are needed for 
each flip in Phase I, and no split is used in Phase II, implying that the total number of splits in 
the algorithm is 0{k). Thus, there are totally 0{n + k) operations on the splittable queues in the 
entire algorithm. By Lemma [3l the total time for performing all k flip operations is 0(n + k). 

In addition, as shown by Lemma [8] (given in Appendix IB.ip . the total time of the preparing 
procedure (used only in Phase I) over the entire algorithm is 0(n + A;). 

We conclude that all k flips can be performed in 0(n + k) time (after the initial good pseudo- 
triangulation To is computed in 0{n + /ilog/i) time). At any moment of the algorithm, the space 
needed is for storing the current good pseudo-triangulation and all splittable queues, which is 0(n). 

If we incorporate the needed graph information into the above algorithm, then the relevant 
visibility graph G can be built in the same amount of time. We thus have the following result. 

Theorem 1 The free bitangents among h pairwise disjoint convex splinegons of totally n vertices 
can be computed in 0{n + k + hlogh) time and 0{n) space, where k is the number of all free 
bitangents. The relevant visibility graph G can also be constructed in 0{n + k + /ilog h) time. 

It remains to prove the key lemma (Lemma (H), which is one major contribution of this paper. 

As in the PV algorithm [30\ I31j . only a constant number of enqueue sequences are involved 
in Phase II for each flip and each enqueue sequence is on either a free bitangent or a boundary 
portion of one single obstacle (see Phase II for more details of this). For the PV algorithm, since 
every obstacle is of 0(1) complexity, each enqueue sequence can be implemented as 0(1) enqueue 
operations. Consequently, Lemma H] easily follows in \30\ I31j. For our problem, however, since the 
complexity of an obstacle can be Q{n), each enqueue sequence may take as many as i}{n) enqueue 
operations. Thus, the simple proof for the PV algorithm apparently does not work for our problem. 
To prove the key lemma, we instead conduct the analysis in a global fashion, as follows. 

As one of our key proof ideas, we define a new concept "reverse" , which was not used in the 
previous analysis [30^ [3T] . Consider a flip on a free bitangent b in Phase I (everything here, such as 
6*, b'j^, R, L, Tji, etc., is defined in the same way as in Sections 13.41 and l3.5p . Let p be a point on 
dT of a pseudo-triangle T in the current good pseudo-triangulation T before the flip such that p 
lies on an obstacle P. Let h (p) be the directed tangent line of T at p. Suppose after the flip of b, p 
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lies on dT' of a new pseudo-triangle T' T); let hip) be the directed tangent line of T' at p. By 
the definition of tangent lines of a pseudo-triangle, both li{p) and /2(p) are tangent to the obstacle 
P at p. Since by our assumption, the boundary of each obstacle is smooth, hip) and hip) lie on 
the same undirected line. But, it is possible that and hip) have opposite directions (e.g., if p 
is on Tr \ {Tailib'fi)} in Case 1; see Fig. [8]). When this occurs, we say that the point p is reversed 
due to the flip of b. Note that p can be reversed only if the pseudo-triangle T is either R or L. 

For example, in Case 1 (resp.. Case 2), all points on Tji except the endpoint Tailih'^ (resp., 
Headib'j^)) are reversed (see Fig. [8]). In Case 3, all points on xrp* \ {p*} (here, xr = Headib)) are 
reversed (note xpip* is part of Tji). Of course, the algorithm does not do the "reverse" explicitly. 

For any atom, if it is (part of) an elementary curve, then we say that it is reversed if all its 
interior points are reversed; if it is a bitangent t, then it is reversed if the direction of t is reversed. 

Our overall proof strategy is to associate the enqueue operations in Phase II of the entire 
algorithm with different "classes" of operations and prove a bound for each such class. For this, 
we denote by n^; the number of enqueue operations in Phase II of the entire algorithm, by ng the 
number of all reversed atoms in the entire algorithm, by the number of dequeue operations 
in Phase I of the entire algorithm, and by ns the number of certain special enqueue operations in 
Phase II of the entire algorithm (which are defined in Appendix IB.3p . Recall that k = \B\. 

Then, proving Lemma H] is to show ue = Oin + k). To this end, we prove that he < tlq + n^) + 
ns + k and ng = 0(n + k), nn = Oin + k), and ns = Oin + k). 

The detailed proof is given in Appendix |Bl which is organized as follows. In Appendix IB.H 
we prove ng = Oin + k). To this end, we prove that any point on any obstacle boundary can be 
reversed at most once in the entire algorithm. We show (in Observation [2|) that the total number 
of all atoms involved in the algorithm is 0(n + A;). We also show (in Lemma [5]) in Appendix 
IB. II that the total running time of the preparing procedure in Phase I of the entire algorithm is 
0(n + A;). In Appendix IB. 21 we prove n^ = Oin + k). To this end, we prove that every atom can 
be dequeued at most 0(1) times in Phase I of the entire algorithm. In Appendix IB.3| we prove 
nE ^ nq + nD + ns + k and ns = 0(n + fc). We show that for each enqueue operation in Phase II, 
say on an atom A, A must belong to one of the following cases: a reversed atom, an atom dequeued 
in Phase I, the current enqueue on A is a special enqueue operation, a free bitangent in B. 

The proof in Appendix [Bl which uses many new observations and analysis ideas, is long and 
technically difficult and complicated. Nevertheless, it does provide lots of insights into the problem 
and explores many essential properties, which may help deal with other related problems as well. 
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Appendix 



A Phase II of a Flip Operation 

In Phase II, our task is to compute Awake[T] and Asleep[T] for each T G {R',L'}. Specifically, 
after Phase II, Awake[T] should represent xrqr of dT and Asleep[T] should represent wtPt of dT, 
where wt = Qt if Qt S Vtzt and wt = yr otherwise. Note that since our definition of Asleep is 
different from that for the PV algorithm [30^ I31j. our algorithmic procedures in Phase II also differ 
from those in the PV algorithm. 

Recall that after Phase I, AwakeMax[R] (resp., AwakeMax[L]) represents p*qR (resp., q*qL) 
in all three cases. In Case 1, AwakeMin[R] (resp., AwakeMin[L]) represents xjiq'^ (resp., XLq'i), 
and in Cases 2 and 3, AwakeMin[R] = AwakeMin[L] = 0. Recall that p* £ xr/Ir and q* S SfiJlL- 

We only show how to compute Awake[R'] and Asleep[R']. The case for L' can be handled 
similarly. We discuss for Cases 1, 2, and 3 individually. No split operation is needed in Phase II. 
Note that after computing b* = ip{b) in Phase I, as in the PV algorithm, the six new cusps, i.e., 
XT, Ut, and zt for T £ {R',L'}, can be determined immediately in 0(1) time. Hence, we assume 
all of them are already known. 

Recall that a splittable queue, e.g., Awake[T] for a pseudo-triangle T, represents a list of 
consecutive atoms on dT. We define the head and tail of the queue such that if moving from the 
head to the tail along the list, it is clockwise on dT. Note that the splittable queue allows to 
enqueue an atom at either the head or the tail of the queue. Recall that pT is the basepoint of T. 
For two points pi and p2 on dT, if pi G prpi, we say pi is before p2; otherwise, pi is after p2 (when 
pi = P2, Pi is both before and after P2). In the following discussion, when a point p is on both dR' 
and dL (resp., dR), we sometimes do not differentiate whether p is on dR' or dL (resp., dR). 



Algorithm 1: Construction of Awake and Asleep of R' for Case 1 
Output: Awake[R'] = xj^qjii and Asleep[R'] = wrPpr' 

1 Awake[R'] AwakeMin[R] ; /* done for Awake[R'] */ 

2 Asleep[R'] ^ AwakeMax[L] ; 

3 Enqueue zlPr/ to the tail of Asleep[R'] ; /* charge the enqueue to Q */ 

4 Enqueue b* to the head of Asleep[R'] ; /* charge the enqueue to B */ 

5 if uri = p* then 

6 Do nothing ; 

7 else/* yp,, = yR */ 

8 Starting at p* , enqueue atoms along dR' counterclockwise to the head of Asleep[R'] until 
either or yR is encountered ; /* charge the enqueue to D */ 

9 end 



In Case 1, recall that bRi = 6^, R' = Rtri{b'^, pRi = Tail{b'^), and qRi = q'^. Note that 
xr' = xr, yR' is yR or p* , and zri is yi or q* . The pseudocode is in Algorithm [H (Note that every 
"enqueue" in the pseudocode is commented by "charge the enqueue to The comments are used 
for analysis in Section [B. 31 and can be ignored at the moment.) The details are discussed below. 

The goal of computing Awake[R'] is to let it represent XR/qRi. Since in Case 1, AwakeMin[R\ = 
xrq'ri and xri = xr and qRi = q'^, by setting Awake[R'] = AwakeMin[R\, we are done. 

The goal of computing Asleep{R'\ is to let it represent wrPpr' ■ Note that in Case 1, AwakeMax[L] 
= q*ZL (since qi = zl) and p* G q'RqR- Thus p* is after (= qRi). Regardless of whether yR' is yR 
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Figure 9: Illustrating Case 2.2 in Phase II: xr = Head{b). 

or p*, p* is after y^/. Thus, p* is after wr', which is either qj^/ or y/j/. Clearly, zl E zjypji'. Hence, 
6* is asleep on dR', and q*ZL, which is stored in AwakeMax[L], is part of Asleep[R']. Thus, we first 
set Asleep[R'] = AwakeMax[L\. Since p^ = Tail{h'^, we enqueue zlRri to the tail of Asleep[R']. 
Now Asleep[R'] represents q*PR>. Since h* is asleep, we enqueue h* to the head of Asleep[R']. Now 
Asleep[R'] represents p*Pr'- Note that y^' is p* or y/j. 

• If i/ij' = p* , then since p* is after = g/j/ , wr' = yji' = p* . Thus we are done with computing 
Asleep[R']. 

• If yji' = yn, then starting at p*, we enqueue atoms along dR' counterclockwise to the head 
of Asleep[R'] until we first encounter or y^. Note that the first point of q'j^ (= qR') or yR 
(= y/j') thus encountered is wr'. 

This completes the construction of Asleep[R']. 



Algorithm 2: Construction of Awake and Asleep of i?' for Case 2.1 
Output: A?i;aA;e[i2'] = XR/qji' and Asleep[R'] = wr>Pr> 

1 As/eepf-R'] ^ ; /* done for Asleep[R'] */ 

2 yl?i;aA;e[i2'] Atua/ceMax [L] ; 

3 Enqueue zlXr of dR to the tail of Awake[R'] ; /* charge the enqueue to Q */ 

4 Enqueue h* to the head of Awake[R'] ; /* charge the enqueue to B */ 

5 if xri = p* then 

6 Do nothing ; 

7 else/* XR' = yR */ 

8 Enqueue yijp* to the head of 74tDaA;e[i?'] ; /* charge enqueue to D */ 

9 end 



In Case 2, recall that bR' = b'j^ and R' = Ltri{b'j^). Note that xr' is yR or p*, yRi is or q* , 
and z/j/ is zl or x/j. Depending on whether xr is Head{b'j^) or Head{b), there are two subcases. 
Case ^.i." = Head{b'j^) (see Fig. ED and Case 2.2: xr = Head{b) (see Fig. [9]). In Case 2.1, 
zr' = XR and in Case 2.2, zri = zl- 

Case 2.1: xr = Headib'^ (see Fig. [8]). In this subcase, Head{b'j^) = pRi = qRi = xr = zri . 
The pseudocode is in Algorithm [2l The details are discussed below. 

The goal of computing Awake[R'] is to let it represent XRiqRi. Recall that AwakeMax[L] 
represents q*ZL on dL. Regardless of whether xr' is yR or p*, both b and the portion q*ZL are 
awake on dR' , i.e., part of XR/qRi . Thus, we first set Awake[R'] = AwakeMax[L] and then enqueue 
zlxr of dR to the tail of Awake[R']. We also enqueue b* to the head of Awake[R']. Now Awake[R'] 
represents p*qR'. Recall that xr' is either yR or p*. If x/j' is p*, then we are done with computing 
Awake[R']. Otherwise {xr' = yR), we enqueue yRP* to the head of Awake[R']. 
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This completes the construction of Awake[R'] for Case 2.1. 

Now we compute Asleep[R'] which represents wj^/pR'. Due to qr/ = pri = zr/, wri is q^t. Hence 
wr^Pr' is only a point, which is not essential to our algorithm. We simply set Asleep[R'] = 0. 
This completes the construction of Asleep[R'] for Case 2.1. 

Algorithm 3: Construction of Awake and Asleep of R' for Case 2.2 
Output: AwakelR'] = XR>qjii and Asleep[R'] = wrPpr' 

1 Determine whether q^i E Pr'P* or p* E PR/Qji' on dR' by checking the position of pri with 
respect to the line containing b* ; 

2 if qri E pr'P* then 

3 Starting at xri, enqueue the atoms to the tail of Awake[R'] along dR' clockwise until qri 
is found ; /* charge enqueue to D */ 

4 else/ * p* e Pr/Qr' */ 

5 Awake[R'] ^ AwakeMax[L] ; 

6 Dequeue atoms from the tail of Awake[R'] until qr' is found ; 

7 Enqueue b* to the head of Awake[R'] ; /* charge the enqueue to B */ 

8 if xr' = p* then 

9 Do nothing ; 

10 else/* xr' = ur */ 

11 Enqueue UrP* to the head of Awake[R'] ; /* charge the enqueue to D */ 

12 end 

13 end 

/* done with Awake[R'] , the following is for Asleep[R'] */ 

14 Asleep[R'] ^ Asleep[L] ; 

15 Enqueue xrPri to the tail of Asleep{R'\ ; /* charge the enqueue to Q */ 

16 if wl = ul then 

17 Do nothing ; 

18 else/* Wl = qi */ 

19 Starting at ql, enqueue atoms along dR' counterclockwise to the head of Asleep[R'] until 
either qRi or yRi is found ; /* charge the enqueue to S */ 

20 end 



Case 2.2: xr = Head{b) (see Fig. [9]). Recall that zri = z^. The pseudocode is in Algorithm [3l 
The details are discussed below. 

We first compute Awake[R'] which represents XRtljRi. Clearly, xri is before qRi on dR' . Since 
we know b* and the basepoint pri of i?', by checking the position of pri with respect to the line 
containing 6*, we can determine whether qRi is before p* on dR' in 0(1) time. Depending on 
whether qRi is before p* on dR' , there are two subcases. 

• qRi is before p* on dR' , i.e., qRi E PR'P* ■ Recall that xri is either yR or p* in Case 2, and xri 
is before qRi on dR' . Since qRi is before p* , xri is before p* , and thus x^j' = yR. 

We first set jdwafcefi?'] = 0, and then starting at xri, we enqueue the atoms to the tail of 
Awake[R'] along dR' clockwise until we find qRf. Now Awake[R'] = XR^qR'. 

• qR' is after p* on dR', i.e., p* E PifqRi. Then clearly, q* E PR^qR'. As can be easily seen, g/j/ is 
before qi on (see Fig. [9]). Hence qRi lies on q*qL of 5L, which is stored in AwakeM ax[L]. 
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b'R 

Figure 10: Illustrating Case 3 in Phase II. 

We set AwakelR'] = AwakeM ax[L], and then dequeue atoms from the tail of Awake[R'] until 
we find qji'. Now Awake[R'] represents q*qR'. We then enqueue b* to the head of Awake[R']. 
Now Awake[R'] represents p*qR'. 

Recall that xr' is either yR or p*. If xri = p*, then we are done with constructing Awake[R']. 
Otherwise {xri = ur), we enqueue vrp* to the head of Awake[R']. 

This completes the construction of Awake[R'] for Case 2.2. 

Next, we compute Asleep[R'] (for Case 2.2), which represents wr/Pr'- Recall that wri is qR' 
or ur', and i/r' is either q* or y^. Below, we will show (implicitly) in a constructive manner that 
wij)L, which is stored in Asleep[L], is part of wr/Pr'. The point wl is either q^ or yi. We need to 
discuss the relations among the positions of q*,yL,QR', and q^. Recall qRi is before qi on dR' (see 
Fig. E]) and q* is before qi on dL. 

We first set Asleep[R'] = Asleep[L], which represents wlPl on dL. Recall that pL = Head{b) = 
xr and PR' = Headih'p^. We then enqueue xrPr' to the tail of Asleep[R']. Now Asleep[R'] 
represents wlPr'- Depending on whether wl is yi or qi, there are two subcases. 

• If Wl = yii then is after q^ on dL. Since q* is before on dL, q* is before yL on 9L, 
implying y/j' = y^. Because yL is after qL, qL is before yR' (= yi). Since qr' is before qL 
on and qL is before y/j' on dR' , qR' is before yRi on Thus, it;/^/ is yR', which is yL 
(= Wl)- Therefore, we are done with the construction of Asleep[R']. 

• If Wl = qL, then yL is before qL on dL and thus Asleep[R'] currently represents qjjpRi- Next, 
starting at qL, we enqueue atoms along dR' counterclockwise to the head of Asleep[R'] until 
we first find qRi or encounter yR' (regardless of whether yRi is yL or g*). Note that the first 
point of qRi and yR' thus encountered is wr'. Now Asleep[R'] represents wr/Pr' and we are 
done with the construction of Asleep[R']. 

This finishes the construction of Asleep[R'] for Case 2.2. 

We thus complete the construction of both Awake[R'] and Asleep[R'] for Case 2. 

In Case 3 (see Fig. [TOl) . = b* , qR' = pR' = p* , and pL = Head{b) = xr. Note that xr' is yL 
or g*, yR' = zl, and z/j' = p*. The pseudocode is in AlgorithmSl The details are discussed below. 

We first compute Awake[R'], which represents XR/qR'. Below, we will show (implicitly) in a 
constructive manner that wlPl, which is stored in Asleep[L], is part of XR/qR' ■ We set Awake[R'] = 
Asleep[L] and then enqueue xrp* (= piJiRi) to the tail of Awake[R']. Now Awake[R'] represents 
wlQr'- Recall that wl is either yL or qL, and qL is after q* on dL. 

• If Wl = yL, then qL is before yL on dL. Thus ?/l is after q* on 5L and x^' is yL (= if^i)- We 
are then done with computing Awake[R']. 
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Algorithm 4: Construction of Awake and Asleep of R' for Case 3 



Output: Awake[R'] = XR^qR' and Asleep[R'] = wb>Pr> 

1 Asleep[R'] ^ ; /* done for Asleep[R'] 

2 Awake[R'] j4s/eep[L] ; 

3 Enqueue xrp* to the tail of Avjake[R'] ; /* charge the enqueue to Q 

4 Enqueue h* to the head of Awake[R'] ; /* charge the enqueue to B 

5 if wl = Vl then 

6 Do nothing ; 

7 else/* Wl = qi 

8 Enqueue xr^Ql to the head of Awake[R'] ; /* charge enqueue to S 

9 end 



*/ 
*/ 



*/ 



*/ 
*/ 



• If Wl = Ql, then ql is after ul on dL, i.e., i/l G xlQl- If a^i?' = 9*1 then g* is on i/lZl', 
otherwise, xr' = ul- In either case, x^i^ hes on yiJlL- We then enqueue xr/Ql to the head 
of Awake[R']. Now 74u)a/ce[i?'] represents xs/qRi and we are done with computing Awake\R'\. 

This completes the construction of Awake[R'] for Case 3. 

We next compute Asleep[R'] (= wrPpr'). Recall that pRi = qui = p* = zri . Thus wri = qRi, 
and wnipR' is only a point, which is not essential to our algorithm. We simply set Asleep[R'] = 0. 
This completes the construction of Asleep[R'] for Case 3. 

B The Proof of the Key Lemma (Lemma [4]) 

This section is devoted to proving the key lemma (Lemma H]) . 

Recall that proving Lemma [His to show ue = 0{n + k). To this end, our strategy is to prove 
that ue < rag + ud + ns + k and nq = 0{n + k), tid = 0{n + k), and ns = 0{n + k). 

The rest of this section is organized as follows. In Section [B.ll we prove nq = 0{n + k). In 
Section IB. 21 we prove ud = 0{n + k). In Section IB. 31 we show ue < nq + no + ns + k and 
ns = 0{n + k). In addition, we show in Section [B.ll that the total running time of the preparing 
procedure in Phase I over the entire algorithm is 0{nq), which is 0{n + k). 

B.l Bounding the Number of Reversed Atoms (i.e., uq = 0{n+ k)) 

Consider a flip operation on a free bitangent b in Phase I. Everything here is defined in the same 
way as in Sections 13.41 and [ 3.5[ e.g., 6*, 6^, R, L, R', L' , Tr, etc. 

In Case 1 (resp.. Case 2), all points on Tr except the endpoint Tail{b'j^) (resp., Head{b'j^)) 
are reversed^Jsee Fig. [8|). In Case 3, all points on xrp* \ {p*} (here, xr = Head{b)) are reversed 
(note that xrp* is part of Tr). Observe that Vr is an obstacle arc. The following observation is 
self-evident. 

Observation 1 After the flip of b, let a be the reversed portion on dT , with T G {R,L}. Let T' 
be R' (resp., L' ) if a lies on dR' (resp., dL'). Then the following properties hold. 

1. a is an obstacle arc. One (resp., the other) endpoint of a is an endpoint of b (resp., bx')- 

2. Let ab (resp., ab^,) be the endpoint of b (resp., bx') on a. From at to a^^, along a, it is 
counterclockwise with respect to the obstacle on which a lies. We call a^, the obstacle-ccw-start 
endpoint of a. 
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3. Every point on a except the endpoint Ob^, is reversed due to the flip of b. For any point 
p £ a \ {ofey, }, a is called the hosting arc of p. 

Note that only points on dR or dL can be reversed due to the flip of b. An important property 
of the reversed portions after every flip is given in the lemma below. 

Lemma 5 After a flip operation on b, suppose a is the reversed portion on dT of a pseudo-triangle 
T, T £ {R, L}. Then no interior point of a can be an endpoint of any bitangent in B. 

Proof: First, by Lemma [21 the direction of any free bitangent t £ B can be reversed only by a 
flip operation on t. Thus, after the flip of b, the direction of b is reversed, whereas the direction of 
any other free bitangent in B does not change. By Observation [H a is an obstacle arc, say, on an 
obstacle P. 

Assume to the contrary that there is an interior point g of a that is an endpoint of a free 
bitangent t £ B. Note that t ^ h since an endpoint of b is an endpoint of a by Observation [TJ 
By the definition of good pseudo-triangulation, the direction of t is compatible with its endpoint q 
before the flip of b. Recall that the direction of t is compatible with its endpoint q if the directed 
tangent line of T at q, denoted by li{q), has the same direction as t. 

Suppose a lies on a pseudo-triangle T' right after the flip of b. Let hiq) be li{q) but with the 
reversed direction. Then both li{q) and hiq) are tangent to the obstacle P at q. Since q is reversed 
due to the flip of 6, hiq) is the directed tangent line of T' at q right after the flip of b. As a free 
bitangent, t's direction does not change after the flip of b t), the direction of hiq) is opposite to 
that of t, making t not compatible with its endpoint q after the flip of b. But this contradicts with 
the definition of good pseudo-triangulation. Thus q cannot be an endpoint of any free bitangent in 
i3, and the lemma follows. □ 

The next lemma is critical. 

Lemma 6 During the topological flip algorithm, any point on the boundary of any obstacle in V 
can be reversed at most once. 

Proof: By Lemma [21 the direction of any free bitangent t £ B can be reversed only by a flip 
operation on t. Since the algorithm does not flip any bitangent more than once, the direction of 
each free bitangent is reversed only once throughout the algorithm. 

Consider a fiip of a free bitangent 6 of a good pseudo-triangulation T. Suppose a point a on 
dT of a pseudo-triangle T is reversed for the first time due to the flip of b. Below we prove that a 
cannot be reversed again. We only discuss the case when a lies on dR (the case on dL is similar). 

Let p = Head{b); let q = Tail{b'^ in Case 1, g = Head{b'^ in Case 2, and q = Tail{b*) in Case 
3, respectively (see Fig. [8]). Due to the flip of 6, in each case, the reversed portion is j3g \ {q}. Also, 
the direction of b is reversed. By Observation [H pq is an obstacle arc, say, on obstacle P. Assume 
to the contrary that later a point a £ pq \ {q} is reversed for the second time. Note that a cannot 
be p since otherwise b would not be compatible with p unless b is reversed twice. 

At the second reversal of a, let p'q' be the hosting arc of a with p' being the obstacle-ccw-start 
endpoint (see Fig. [TT]) . By Observation [H p'q' is an obstacle arc, say, on obstacle P' . Since the 
point a is on both pq £ P and p'q' £ P', a lies on both P and P' . Therefore, P = P' since our 
obstacles in V are pairwise disjoint. Also by Observation [H both p' and q' are endpoints of some 
free bitangents in B. 

Since p' is the obstacle-ccw-start endpoint of p'q', by Observation [H when moving from p' to q' 
on p'q', it is counterclockwise with respect to P. Similarly, when moving from p to q on pq, it is 
counterclockwise with respect to P. Since a is on both pq and p'q', there are three possible cases: 
(i) p is an interior point of p'q', (ii) p' is an interior point of pq, and (iii) p = p'. We argue below 
that any case can not occur. Consequently, the point a cannot be reversed again. 
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Figure 11: Illustrating the proof of Lemma (6] 

As the point p is an endpoint of b, by Lemma O p cannot be an interior point of p'q'. Similarly, 
p' cannot be an interior point of pq. For case (iii), since the point p' is the obstacle-ccw-start 
endpoint of p'q', p' is reversed. If case (iii) occurs, p' (= p) is the endpoint of 6, and thus the 
direction of b has to be reversed again since otherwise b would not be compatible with the reversed 
p'. But this contradicts with the fact that b cannot be flipped twice in the algorithm. 

The lemma thus follows. □ 
The following lemma bounds the number of all different atoms involved in the algorithm. 

Observation 2 The total number of different atoms involved in the algorithm is 0{n + k). 

Proof: In our problem, an atom is of one of the following three types: (1) A free bitangent, (2) 
an elementary curve, and (3) a portion of an elementary curve. It is easy to see that the number 
of type (1) atoms is 0{k) and the number of type (2) atoms is 0{n). To prove this observation, it 
suffices to show that the number of type (3) atoms is 0{n + k). 

After initialization, in the initial good pseudo-triangulation To, there are 0{n) type (3) atoms. 
Subsequently in the algorithm, each new type (3) atom is produced only due to a flip operation. 
Note that after every flip operation, at most 0(1) new type (3) atoms can be produced. Since there 
are k flips, plus those in To, the number of type (3) atoms is bounded by 0(n + k). □ 

We then have the following result. 

Lemma 7 The number of reversed atoms in the entire algorithm (i.e., uq) is 0{n + k). 

Proof: By Lemma [6l each atom can be reversed at most once. Consequently, Observation [2] leads 
to the lemma. □ 

Lemma 8 The overall running time of the preparing procedure in Phase I of the entire algorithm 
is 0{n + k). 

Proof: Recall that the preparing procedure is to determine whether Case 1 occurs and find b'^ if 
Case 1 holds. Since the three cusps of R are maintained by the algorithm, xr is known. We walk 
from Head{b) along dR clockwise until first encountering either b'j^ or xr. If we encounter b'j^ first, 
then it is Case 1 and we have 5^ as well. Otherwise, it is either Case 2 or Case 3. 

It is easy to see that the running time of the above walking is proportional to the number 
of atoms of the reversed portion due to the fiip of b. Therefore, the total time of the preparing 
procedure for the entire algorithm is at most proportional to the total number of reversed atoms 
in the entire algorithm, i.e., Oinq), which is 0{n + k) by Lemma [71 The lemma thus follows. □ 

B.2 Bounding the Number of Dequeue Operations in Phase I (i.e., no = 0{n+k)) 

Our goal in this section is to prove Lemma [9] below. 
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Lemma 9 In Phase I of the entire algorithm, the total number of dequeue operations (i.e., no) is 
0{n + k). 

To prove Lemma [9l we will first prove the following lemma (Lemma llOp . which states that any 
atom can be dequeued in Phase I of the entire algorithm at most twice before its (only) reversal 
and can be dequeued in Phase I at most twice after its reversal (if any). Note that Lemma [TOl does 
not include the dequeue operations in Phase IL 

Lemma 10 For any atom A in the algorithm, regardless of whether it has been reversed previously, 
A can be dequeued in Phase I at most twice before its next reversal ( if any ) in the algorithm. 

By Lemma [6l every point on &P can be reversed at most once. This implies that, by Lemma 
[TUt every atom can be dequeued in Phase I at most four times in the entire algorithm. Since the 
total number of atoms in the algorithm is 0{n + /c), Lemma [9] follows. 

The rest of this section gives the proof for Lemma [TOl The proof, which is quite long and 
technically complicated, is based on many new geometric observations and analysis techniques. 

The following lemma (proved in \?>Q\ I31j ) will be repeatedly referred to by our analysis later. 

Lemma 11 For any pseudo-triangle T of a good pseudo-triangulation, (1) if yx lies on 

xtQt of dT , then yrOr is an obstacle arc, and (2) if zt i^PT, then zxpr is an obstacle arc. 

For any bitangent t G B{T) of a good pseudo-triangulation T in the algorithm, we view t as 
defining two atoms, one for Rtri{t) and the other for Ltri(t), and we say that these two atoms 
are defined by t. Similarly, we let every point on t have two copies that belong to Rtri{t) and 
Ltri{t), respectively. Thus, for any point p G dT for a pseudo-triangle T in T, if p is on a bitangent 
t G B{dT), then p refers to the copy of the point on t that belongs to T. In this way, every point 
on dT belongs to exactly one pseudo-triangle in T, i.e., the pseudo-triangle T. Since there are k 
bitangents in B, by Observation [2l the total number of atoms is still 0{n + k). When an atom A 
is (part of) an elementary curve, we also say A is an obstacle arc. 

Before presenting the main proof, we give some observations, which will be useful later. 

B.2.1 Some Observations 

Observation 3 Consider a pseudo-triangle T in a good pseudo-triangulation at a (time) moment 
^1 of the algorithm. For any point p E dT, let li (p) be the directed tangent line of T at p. Suppose 
later at a moment ^2 of the algorithm, the point p lies on dT' for a pseudo-triangle T' , and let 
hip) be the directed tangent line of T' at p. Then, li{p) and hip) He on the same undirected line. 
Further, if p is not reversed during the time period from ,^1 to ^2, then li{p) has the same direction 
as hip); otherwise, hip) and hip) have opposite directions 

Proof: The observation can be easily proved by the definition of a directed tangent line of a 
pseudo-triangle. If A is defined by a free bitangent, then hip) and hi?) both lie on the same 
undirected line that contains A. If A is an obstacle arc, then hip) and hip) also both lie on the 
same undirected line that is tangent to the obstacle where p lies. Thus, in any case, hip) and hip) 
lie on the same undirected line. 

If the atom A has not been reversed during the time period from ^1 to .^2) then clearly hip) and 
hip) are of the same direction. Otherwise, by Lemma [6l p is reversed only once during the time 
from ^1 to ^2, and therefore, hip) and hip) have opposite directions. □ 

Observation 4 Consider a pseudo-triangle T in a good pseudo-triangulation at a (time) moment 
^1 of the algorithm. For any two points p and q on dT, let lip) and /(g) be their corresponding 
directed tangent lines of T. Then the following properties hold. 
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1. If the l{q)-slope ofl{p) is less than vr, then the l{q)-slope of l{p) is no bigger than the pt-slope 
of l{p) in T. 

2. If at a later moment ^2 of the algorithm, the point p lies on dT' of a pseudo-triangle T' 
(T' = T is possible) such that the same directed l{p) is the directed tangent line of T' at p 
then the pt-slope ofl{p) at the moment ^1 is no smaller than the pt-slope ofl{p) at the moment 

6- 

Proof: For the first part, observe that the pt-slopes of ah points on dT {l{p) and l{q) included) 
are no more than vr. Suppose the Z(g)-slope of l{p) is less than vr. If we rotate bx around its 
tail counterclockwise, we will encounter first the direction of l{q) and then the direction of l{p) 
(otherwise, the pt-slope of l{q) would be larger than vr). This means that the pt-slope of l{p) is the 
sum of the pt-slope of l{q) and the /((?)-slope of l{p). The first part thus follows. 

The second part can be proved by a simple induction. Recall that if the point p lies on a 
bitangent in B{dT), then p refers to the copy of the original point that belongs to T. If T' = T, 
then it is easy to see that the property holds. Now consider the first flip operation after which the 
point p lies on dT' of a pseudo-triangle T' ^T such that l{p) is still the directed tangent line of T' 
at p. Since T' 7^ T, the above flip must be on the free bitangent 6t- 

Let ^'i (resp., ^2) be the moment right before (resp., after) the flip of hx- At the moment 
one endpoint of 6^' must be on dT although bx' rnay be (p{bT)- Recall that due to the flip of bx, 
every free bitangent in ;S \ {6t} does not change direction. So bx' does not change direction due 
to the flip of bx- Consequently, the 6r-slope of bx' on dT at the moment ^[ is the pt-slope of the 
endpoint of bx' that is on dT, which must be less than vr. Note that the fe^-slope of l{p) is the 
pt-slope of l{p) on dT at the moment ^J, which must be less than vr. Similarly, the 67^/-slope of 
l{p) is the pt-slope of l{p) on dT' at the moment which must be less than vr. Therefore, if we 
rotate bx around its tail counterclockwise, we will encounter flrst the direction of bx' and then the 
direction of l{p) (otherwise, the fe^-slope of bx' would be larger than vr at the moment ^^). This 
implies that the iy-slope of l{p) at the moment is no smaller than the ^T'-slope of l{p) at the 
moment i-^-, the pt-slope of l{p) at the moment is no smaller than the pt-slope of l{p) at the 
moment ^2- 

Consider the flrst flip after which the point p lies on dT" of a pseudo-triangle T" 7^ T' such 
that l{p) is still the directed tangent line of T" at p. Let ^3 be the moment right after this flip. By 
a similar argument, the pt-slope of l{p) at the moment ^2 is no smaller than the pt-slope of l{p) at 
the moment ^3. Inductively, the second part holds. □ 

Note that Observation 111^2) above actually tells us that for any point p G dT of a pseudo- 
triangle T in the algorithm, regardless of whether p has been reversed previously, the pt-slope of p 
is monotonically decreasing during the rest of the algorithm until it is possibly reversed. 

B.2.2 The Main Proof of Lemma [H 

In the following proof of Lemma [T0| all dequeue operations refer to those in Phase I. 

Consider the dequeued atoms on dR due to the flip of b for the pseudo-triangle R = Rtri{b). 
In Case 1, the dequeued atoms are on q'^p*; in Cases 2 and 3, the dequeued atoms are on xnp* . In 
all three cases, the dequeued atoms of dR are also on dR' . Recall R' = Rtri{b*) and L' = Ltri{b*). 
The situation on the pseudo-triangle L = Ltri[b) is similar. 

Let A be an arbitrary atom of the good pseudo-triangulation T right before the flip of b. Note 
that A may have been reversed before (and thus cannot be reversed again), but this is not important 
to our proof. Observe that due to the flip of 6, A can be dequeued only if (a) A is awake before 
the flip of 6, and (b) A lies on either dR or dL. Note that b = ba = bL since b is minimal. In other 
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words, suppose A lies on dT of a pseudo-triangle T; if A is dequeued due to the flip of b, then (a) 
T is either R or L, (b) A must be awake on dT (before the flip of b), and (c) b = bx- 

Our main idea for the proof is to show that, if an atom A has been dequeued, then before 
its next reversal (if any), A cannot be awake (and thus cannot be dequeued again). To this end, 
one observation used in the proof is that: For a pseudo-triangle T with T = Rtri{bT) (resp., 
T = Ltri{bT)), if a point p is awake on dT, then the forward (resp., backward) T-view of p has a 
smaller pt-slope than p in T (see Fig. [7]). The proof, however, consists of a lengthy and complicated 
case analysis with considerable details. 

To prove Lemma [TOl we first prove the following statement, which we call Sub-lemma 110(a). 

Sub- lemma 1101 (a). Suppose A is an atom on dT of a pseudo-triangle T with T = Rtri{bT) in a 
good pseudo-triangulation T and A is dequeued due to a flip operation on br- Also, suppose at any 
later moment before the reversal of A, A lies on dT' of a pseudo-triangle T' with T' = Rtri{bj'r) in 
another good pseudo-triangulation T'. Then A cannot be dequeued due to the flip of bx'- 

To prove Sub-lemma llOf a) below, we analyze the three main cases (Cases 1, 2, and 3). 

The proof of Case 1. In Case 1, all dequeued atoms lie on q'j^p*. Let A be an arbitrary atom on 

q'jlP*- Let .^1 be the moment right after the flip of b. Hence the atom A is on dR' at the moment 
^1. Suppose at a later moment ^2 of the algorithm, the atom A lies on dT' of a pseudo-triangle 
T' of a good pseudo-triangulation T' with T' = Rtri{bT') and A has not been reversed since the 
moment ^1. As discussed above, the atom A on dT' can be dequeued only due to the flip of 6t'- 
Thus, during the time between the moment ^2 and the flip of bx', A cannot be dequeued on dT'. 
Note that once it is formed, the pseudo-triangle T' remains unchanged in the algorithm (and thus, 
A is not reversed) until bx' is flipped. Without loss of generality, we let ^2 be the moment right 
before the flip of 6^/. We prove below that A cannot be awake at the moment ^2 and thus cannot 
be dequeued due to the flip of ftr' • In the following discussion, for simplicity, when we mention R' 
(resp., T'), we always refer to the moment ^1 (resp., ^2) unless otherwise stated. 
Let p be an arbitrary interior point on A, i.e., p is not an endpoint of A. 

Let hip) be the directed tangent line of T' at p (at the moment ^2), and li{p) be the directed 
tangent line of R' at p (at the moment ^1). Since p is not reversed during the time period from .^1 
to .^2 5 by Observation [3l li{p) and hip) are the same. Below, we simply use l{p) to refer to both 
h{p) and hip)- Since p is not reversed due to the flip of b, lip) is also the directed tangent line of 
R at p before the flip of b. 

Below, we prove that at the moment ^2) the point p is not awake on dT' (i.e., p does not lie on 
Awake[T'] = xt^Qt') and thus the atom A cannot be awake. Let q G dT' be p's forward T'-view 
point along lip) (at the moment ^2)- Let /(g) be the directed tangent line of T' at q. 

Assume to the contrary that the point p is awake on dT' (at the moment ^2)- Then it immedi- 
ately implies that p lies on AwakelT'] = xjyqx' and the pt-slope of l{p) is larger than that of /(g) 
in T' (see Fig. [7]). Thus, /(g) must cross lip) from left to right. 

Note that right before the flip of b, we have Pl = qi = zi = Headib), i.e., pL = zl is the 
basepoint of L = Ltriib) and qL is the point on dL whose backward L-view (on dL) is pL and 
qL = zl in this case (see Fig. [8|). Let c = Tailib'j^). 

Recall that the point p is an interior point of an atom A on q'j^p*. Recall that is the point 
on dR whose forward R-view is Tailib'j^. If g^ is on a bitangent t G BidR), since p is the interior 
point of t, every point on t can be viewed as g^; in this case, we let g^ be the endpoint of t such that 

t lies entirely on xnq'^. This step can be done when we conduct the split operation on Awake[R\ 
at q'^ (i.e., change the criterion when searching q'^ in Awake[R\). Note that the above requirement 
for q'j^ does not change the running time of Lemma [3l In this way, / (p) must intersect either yLzZ 
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Figure 12: Illustrating the case that the point g is a Figure 13: Illustrating the case that g lies on /(p) before 
(here a £ ulzl)- a, (here a £ i/lzl)- 

on dL or z^c on dR (see Fig. [12] or Fig. [T3|) . 

Let a be the intersection point between l{p) and vlzl or zlc (see Fig. [12] or Fig. [T3]l . Since zlc 
is the reversed portion on dR due to the flip of b, by Observation [T] zlc is an obstacle arc. Since 
ql = zl, by Lemma [TT] the part (= Vlzl) is an obstacle arc. Hence, the point a must be on 

an obstacle, say P. Consider the position of g, which is p's forward T'-view point on dT' along l{p) 
(at the moment ^2)- Since a lies on an obstacle P, q can be either at a or on l{p) between p and 
a (but before a). In the following, we show that either case cannot occur, and consequently our 
assumption that the point p is awake on dT' is not correct. 

If q is at the point a (see Fig. [T2|) . then l{q) is the directed tangent line of T' at a {= q)- Since 
a lies on the obstacle P, l{q) is tangent to P at a. Note that the point a lies on dR' at the moment 
^1 (right after the flip of h). Let li{q) be the directed tangent line of R' at a (at the moment ^1). 
Since a lies on P, li(g) is also tangent to P at a, and therefore, li{q) and ^(g) both lie on the same 
undirected line. There are two subcases to consider: a lies on 22c or on yizZ \ {^l}- Below, we 
show that neither case can occur. 

(i) If a (= q) lies on zlc (a can be zl), since is an endpoint of an atom on q'^p* and p is an 
interior point of the atom A, we have p q'p^ and a ^ c. Since 'zlc is reversed due to the flip 
of 5, by Lemma [6] 'zic will not be reversed again after the moment .^1. Thus, l{q) has the 
same direction as li{q). Recall that l{q) crosses l{p) from left to right. Thus, li{q) crosses 
l{p) from left to right. However, at the moment ^1, both p and q (= a) are on dR' and it is 
easy to see that the pt-slope of l{p) is less than the pt-slope of li{q) in R' . Thus, l{p) must 
cross li{q) from left to right, or equivalently, li{q) must cross l{p) from right to left. Hence, 
we obtain that li{q) crosses l{p) both from left to right and from right to left, which is a 
contradiction. Thus, a cannot lie on z^c. 

(ii) For the subcase when a (= q) lies on yLzZ \ {zl} (see Fig. [T2]) . if a has not been reversed since 
^1, then at the moment the direction of l{q) is the same as /i(g)'s. By a similar argument 
as for the former subcase (i), we can show that a contradiction also occurs. 

In the following, we assume that a is reversed (only once) during the time from ,^1 to ^2- 
Thus, li{q) and l{q) have opposite directions. Since l{q) crosses l{p) from left to right, li{q) 
crosses l{p) from right to left. Since l{q) crosses l{p) from left to right, the /(g)-slope of l{jp) 
must be less than vr. Since both l{p) and l{q) are directed tangent lines of T', by Observation 
m^l), the pt-slope of l{p) in T' is no smaller than the /(g)-slope of l{p) (at the moment ^2)- 

Below, we consider l{q) as a physically directed line that is not associate with any time 
moment. We claim the /(q)-slope of l{p) is larger than the 6i?'-slope of l{p) at the moment .^i. 
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Indeed, the 6i{/-slope of li{q) is the pt-slope of the point q (= a) on dR', which is larger than 
zero and less than vr. Since li{q) and l{q) have opposite directions, the Z(g)-slope of bji/ is less 
than vr and larger than zero. Note that l{p) is the directed tangent line of R'. So the bpt'-slope 
of l{p) is the pt-slope of l{p) in R' (at the moment ^i), which is less than vr. To summarize, 
we have: (1) the /((7)-slope of is less than vr and larger than zero, (2) the bRi-slope of l{p) 
is less than vr, and (3) the l{q)-slope of l{p) is less than vr. Therefore, the /(g)-slope of l{p) is 
the sum of the l{q)-slope of bji' and the i^z-slope of l{p). Since the l{q)-slope of 6^^/ is larger 
than zero, the claim is true. Since the 6^'-slope of l{p) is the pt-slope of l{p) in R' , we obtain 
that the /(g)-slope of l{p) is larger than the pt-slope of l{p) in R' at the moment ^i. 

To summarize what have been deduced above, we have: (i) the Z(g)-slope of l{p) is larger 
than the pt-slope of l{p) in R' at the moment ^i, and (ii) at the moment the pt-slope of 
l{p) in T' is no smaller than the Z((;)-slope of l{p). These imply that the pt-slope of l{p) in R' 
at the moment .^i is smaller than the pt-slope of l{p) in T' at the moment ^2- However, l{p) 
is the directed tangent line of both the pseudo-triangles R' and T'; by Observation |4l(2), the 
pt-slope of l{p) in R' at the moment must be no smaller than the pt-slope of l{p) in T' at 
the moment ^2) which incurs a contradiction. 

Hence, we conclude that a cannot lie on yzZL \ {zl}- 
The above analysis shows that q cannot be at the point a. 

We then discuss the case when q lies on l{p) between p and a (but before a). In other words, 
l{p) intersects l{q) (at q) before a at the moment .^2) as shown in Fig. [T3j Since the interior of the 
line segment pa connecting a and p (lying on l{p)) does not intersect any obstacle, it follows that q 
lies on a directed free bitangent ^2 in B{dT') at the moment ^2- Clearly, t2 has the same direction 
as l{q) and thus t2 crosses l{p) from left to right. We let t2 be a physical copy of t2 (i.e., they are at 
the same location with the same direction) but t2 is not associated with any time moment. Then 
t2 crosses l{p) from left to right as well. Note that the interior of R' is free of obstacles and the 
segment pa is contained in R' . Because pa intersects t2 before a, we claim that there must be a 
(directed) bitangent t'l G B{dR') on xj^p or pzjii such that t2 crosses t'l from left to right. This 
claim is proved in the next paragraph. 

Recall that dR' consists of three convex chains, i.e., xr^Jr/, yifzjii, zr^xr'. Note that XRppVJ 
0R' = ^i^r/iJr' U Ur'Zr'- Since the interior of R' is free of obstacles and pa (which is contained in 
R') intersects the directed free bitangent t2 (at q) before a, t2 must cross dR' somewhere, at xr/p 
or pzR/ (and possibly at other locations of dR'). We discuss below the subcase when t2 crosses 
xri~p (the other subcase can be analyzed similarly). Let w be the first point of xr/~p encountered as 
walking on t2 from q in the direction of t2- Then since t2 is a free bitangent, the point w must lie on 
another free bitangent on XRpp. Below, we show that t2 crosses t.^ from left to right. Note that 
in Case 1, the portion wp of dR' does not contain the basepoint pri (= Tail{b'j^)) of R'. Further, 
wp is to the right of both t2 and l{p). Let w' be the endpoint of the bitangent t^, that lies on wp. 
As a portion of u^, w'p does not contain the basepoint pr' of R' and w'p is to the right of both t2 
and l{p). Suppose we move a point w" from p along w'p to w'; let p{w") be the ray originating at 
w" and shooting in the direction of the directed tangent line of R' at w" (i.e., p{w") is the directed 
half-line of the directed tangent line of R' at w"). Then since w'p does not contain the basepoint 
PR' of R', the direction of p{w") changes continuously as we walk along w'p. In particular, when 
w" is at p, p{w") lies on l{p) and has the same direction as l{p)] when w" arrives at w' , p{w") 
contains and has the same direction as t^. Note that l{p) crosses t2 from right to left. Since the 
direction of p{w") changes continuously for w" G w'p and w'p is to the right of both t2 and l{p), 
during the movement of w" from p to w' on w'p^ the ray piw") always crosses t2 from right to left. 
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In particular, when w" arrives at w' , p{w') crosses t2 from right to left. Since p{w') has the same 
direction as tw and t2 crosses (at w), the directed bitangent tyj crosses t2 from right to left, or 
equivalently, t2 crosses from left to right. Letting t'^ = t^, the claim holds. 

Let t'2 be the version of the bitangent t'^ at the moment ^2 (i-e., t'l and t'2 are defined by the 
same undirected free bitangent but may have different directions). So t'2 has opposite direction to 
t'^ if and only if t'^ is flipped during the time period from (^1 to ^2- Since ^2 crosses i'^, t2 crosses at 
the moment ^2- Recall that at the moment ^2, ^2 is in B{dT'). Thus, since t2 crosses t^, cannot 
be in B{T') for the good pseudo-triangulation T' at the moment ^2- Because t'^ G B{dR') at the 
moment .^1, there must be one and only one flip operation on t'^ during the time from ^1 to ^2) which 
reverses the direction of t'l- Hence, t'^ and t'2 have opposite directions. Since t2 crosses t'l from left 
to right (at the moment ^1) and t2 has the same direction as t2, ^2 crosses t'2 from right to left at 
the moment ^2- However, at the moment ^2, we have ^2 G B[T') and t'2 B{T') for the current 
good pseudo-triangulation T'; by the third property of the definition of good pseudo-triangulation, 
t2 should cross t'2 from left to right. This incurs a contradiction. 

Consequently, the case when q lies on l{p) between p and a (but before a) cannot occur. 

Therefore, our assumption that the point p is awake on dT' at the moment ^2 is not correct. In 
other words, the point p cannot be awake at the moment ^2- Consequently, the atom A cannot be 
awake at the moment ^2 (i-e., right before the flip of 6^/). 

As a summary for Case 1, we conclude that when T' = Rtri{bT'), the atom A cannot be 
dequeued due to the flip of 6^'- This flnishes the proof of Case 1 for Sub-lemma [TUlfa) . 

In Case 2, all dequeued atoms lie on xrp*. Depending on whether xr is Head{h'^) or Head{b), 
there are two subcases. Case 2.1: xr = Head{h'^) (see Fig. [8|) and Case 2.2: Head{b) = xr (see 
Fig. [T^ . For convenience. Case 2.2 will be analyzed after Case 3. 

The proof of Case 2.1. In Case 2.1, let A be an arbitrary atom on xjip* (dequeued due to the 
flip of b), and ^1 be the moment right after the flip of b. Suppose at a later moment ^2 of the 
algorithm, the atom A lies on dT' of a pseudo-triangle T' of a good pseudo-triangulation T' with 
T' = Rtri(bT') and A has not been reversed since the moment ^1. Without loss of generality, let ^2 
be the moment right before the flip of bx'- By a similar analysis as for Case 1, we can prove that 
A cannot be awake at the moment ^2 and thus cannot be dequeued due to the flip of bx'- Below 
we sketch the similarity and (minor) difference between the analysis for Case 2.1 and Case 1. 

Let p be an arbitrary interior point on A, and l{p) be the directed tangent line of T' at p at 
the moment ^2- Note that in Case 2.1, Head{b) = zl = ql = pi still holds (see Fig. [8|). Further, 
although Case 2.1 does not involve with q'j^, the same critical structure for this case as for Case 1 
is that l{p) must intersect either iJlzl on dL or zlxr on OR, both lying on the same obstacle, say 
P. To see this, first, since zlxr is the reversed portion on dR due to the fiip of 6, by Observation 
El zlxr is an obstacle arc. Second, due to qi = zl, by Lemma [TTl the portion yiOL (= Vlzl) is 
an obstacle arc. Let a be the intersection of l{p) with yLzZ or zixr on dP. As in Case 1, the fact 
still holds that p's forward T'-view point q on dT' at the moment ^2 is either at the point a or on 
l{p) between p and a (but before a). Hence, the rest of the analysis simply follows as in Case 1. 

We conclude that in Case 2.1, when T' = Rtri{bT'), the atom A cannot be awake at the moment 
^2 and consequently cannot be dequeued due to the fiip of b^'- This finishes the proof of Case 2.1. 

The proof of Case 3. For Case 3, all dequeued atoms lie on xjip* of OR, which are immediately 
reversed after the flip of b (see Fig. [8|). That is, for each atom A on x^p*, after it is dequeued, it 
is reversed immediately as well. Hence in this case, for each dequeued atom A, it is obviously true 
that A is not dequeued again before its forthcoming reversal (due to the flip of b). 

The proof of Case 2.2. For Case 2.2 (see Fig. [Ml) , all dequeued atoms lie on x^p*. Let c = 
Head{b'^). Note that the portion xrc is reversed due to the flip of b. Similarly to the analysis for 
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Figure 14: Illustrating Case 2.2 in Phase I: xr = Head{b). 



Case 3, it is clearly true that any dequeued atom on xrc is not dequeued again before its reversal. 

For the portion cp* of x^ip*, let A be an arbitrary atom on cp*, and be the moment right 
after the flip of b. Suppose at a later moment ^2 of the algorithm, the atom A lies on dT' of 
a pseudo-triangle T' of a good pseudo-triangulation T' with T' = RtriQ}^') and A has not been 
reversed since the moment ^1. Further, let ^2 be the moment right before the flip of b^'- Next, we 
prove that A cannot be awake at the moment ^2 and thus cannot be dequeued due to the flip of 
bx'- Parts of the proof use similar analysis techniques as for Case 1, which we will only sketch. 

Let p be an arbitrary interior point on A, and l{p) be the directed tangent line of T' at p at the 
moment ^2- For simplicity, we always associate R' (resp, T') with the moment ^1 (resp., (,2)- As in 
Case 1, l{p) is also the directed tangent line of R' at p. 

Let q £ dT' be p's forward T'-view point on l{p), and l{q) be the directed tangent line of T' at q 
(at the moment ^2)- Recall that a point w on dT' is awake if and only if w £ xr/qx' (= Awake[T']). 

Assume to the contrary that the atom A on dT' is awake at the moment ^2 • Then it immediately 
implies that the point p £ A lies on Awake[T'] = xjyQT' and the pt-slope of l{p) is larger than that 
of l{q) in T' (see Fig. [Tj). Thus, l{q) must cross l{p) from left to right. 

We first briefly explain why the argument in Case 1 does not work in this case. Note that l{p) 
must intersect yiyzji' on dR' (see Fig. [H]); let a be the intersection point of l{p) and Vb^zri . In 
Case 1, the corresponding point a must lie on an obstacle and thus p's forward T'-view point q on 
dT' at the moment ^2 is either at the point a or on l{p) between p and a (but before a). However, in 
Case 2.2, the point a may not lie on an obstacle. Consequently, the point q may lie on l{p) beyond 
the point a, i.e., "behind" the chain yifz]^'. This makes the argument in Case 1 not applicable to 
Case 2.2. A new analysis approach is given below. 

Clearly, the point q can be on either an obstacle or a bitangent in B{dT'). In the following, we 
show that neither of these two cases can occur (the analysis techniques are somewhat similar), and 
consequently our assumption that the atom A is awake on dT' is not correct. 

We first consider the case when q lies on a bitangent in B{dT'), denote by tq (at the moment 
^2)- Since tg has the same direction as l{q), tq crosses l{p) from left to right. Let ig(^i) be the 
version of the bitangent tq at the moment ^1 (i.e., tg('^i) and tq are defined by the same undirected 
free bitangent but may have different directions). So tq has opposite direction to t'^{£^i) if and 
only if t'q{^i) is flipped during the time period from ^1 and ^2- Recall that by our general position 
assumption, no three obstacles have a common tangent line. Note that l{p) must intersect Vb/'zrj 
on dR' (see Fig. [14]); let a be the intersection point of l{p) and ynjzji' ■ Depending on the relations 
between the free bitangent tq(^i) and the chain yifzjii, there are further four subcases to consider: 
(i) t'q{^i) is part of yifzR'] (ii) tg(^i) crosses y^^^^R' (i.e., tg(^i) crosses a bitangent in yifz^i since 
ig(^i) is a free bitangent); (iii)^^^^i) does not intersect yjfzji' and q lies on l{p) between p and a; 
(iv) ig(Ci) does not intersect vrizri and q lies on l{p) beyond a. Note that in subcase (ii) above, 
q can be either between p and a or beyond a on l{p)- We prove below that none of these four 
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subcases can occur. 

(i) tg(^i) is part of yjfz]^'. So q lies on yRiz^i of dR' in this subcase. Clearly, the point p lies 
on cy^i of dR' (recall c = Head{b'^)). Recall that = h'^ (in Case 2). Thus, the pt-slope 
of l{p) in R' is less than that of t'^Hi)^ and thus l{p) must cross tg(^i) from left to right. 
Since i'q(^i) G S(5i?') and G B{dT') C 5(r') for the good pseudo-triangulation V at 
the moment .^2, the bitangent t'q{^i) is not flipped during the time period from ^1 to ^2 since 
otherwise tq would not be in B{T'). This implies that tq(Ci) and tq have the same direction. 
Thus l{p) must cross tq from left to right, or equivalently, tq must cross l{p) from right to left, 
contradicting with that tq crosses l{p) from left to right. Hence, this subcase cannot occur. 

(ii) tg(^i) crosses Vb/z^' (i.e., iq(^i) crosses a bitangent on yifzji'). As in Case 1, since the interior 
of R' is free of obstacles, it is easy to see that tq(^i) must also cross zj^fxRi or xj^/tjR' of dR' . 
Note that due to the flip of b, the portion xrc on dR is reversed. By Observation [H xrc is 
an obstacle arc. Note that = zri in Case 2.2 (see Fig. Since pL = Head{b) = xr, by 
Lemma [TTl zijpl = Sr^^Tr is an obstacle arc. Thus, ig(Ci) cannot cross the subchain zj^c of 
zr'Xh' on We let tq be a physical copy of tq (i.e., they are at the same location with 
the same direction) but tq is not associated with any time moment. Then tq crosses cxr' or 
XRiyfii of dR' . By a similar analysis as Case 1, we can show that there must be a directed 
free bitangent t'l G B{dR') on cp or py/j' such that tq crosses t'l from left to right. Again, 
as the analysis for Case 1, this will incur a contradiction with the third property of good 
pseudo-triangulation (at the moment ^2)- Therefore, this subcase cannot occur. 

(iii) tq(^i) does not intersect ytfzRi and q lies on l(j}) between p and a. In this subcase, since the 
interior of R' is free of obstacles, tg(Ci) niust cross zj^iXri or xjyyji' of dR'. The rest of the 
analysis follows that of subcase (ii) above. We conclude that this subcase cannot occur. 

(iv) tg(^i) does not intersect ytfz^' and q lies on l{p) beyond a. Note that g / a since tq(^i) does 
not intersect yni'zR' ■ Clearly, a lies on a free bitangent in yj^/^' ; let t'^ denote this bitangent 
(at the moment ^1). Note that l{p) crosses t'^ from left to right since the pt-slope of l{p) 
is less than that of t'^ in R'. Let ta{(,2) be the version of at the moment ^2 (i-e., t'^ and 
ta{^2) are defined by the same undirected free bitangent but may have different directions). 
Since q ^ a, l{p) intersects ia(^2) (at a) in the interior of the pseudo-triangle T', and thus 
ta{^2) cannot be in B{T') for the good pseudo-triangulation T' at the moment ^2- Since 
t'g^ E B{dR'), there is one and only one flip on t'^ during the time from ^1 to ^2- Thus, t'^ 
and ta{S,2) have opposite directions. This implies that l{p) crosses ta{S,2) from right to left, or 
equivalently, ta{S,2) crosses l{p) from left to right. Consider the two endpoints of ia(^2) with 
respect to dT'. There are three possibilities: (a) Both endpoints of ia(^2) are on dT'; (b) 
only one endpoint of ia(?2) is on dT'; (c) neither endpoint of ta(?2) is on dT'. We show below 
that none of the above possibilities can occur, and thus this subcase cannot occur. 

(a) Both endpoints of ia(^2) are on dT'. Note that for any free bitangent t* with t* 
B{dT'), the two endpoints of t* cannot be both on dT' . Since ta{S,2) B{T'), we have 
taiS,2) B{dT') and thus this possibility cannot occur. A detailed proof is given below 
by contradiction. 

Assume to the contrary that both endpoints of ta{(,2) are on dT' . Let a' and a" be 
the two endpoints of ta{(,2) such that the interior of a' a" on dT' does not contain the 
basepoint px' of T' . Denote by P' the obstacle where a' lies. Let / be the directed line 
containing ia(^2) with the same direction as ta{£,2)- Then I is tangent to P' at a'. Since 
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a' is on dT' , by the definition of tangent lines of a pseudo-triangle, I is also tangent to T' 
at a' . Because I has the same direction as ta(C2)) I is the directed tangent line of T' at 
a' (at the moment ^2)- Similarly, I is also the directed tangent line of T' at a" . Hence, 
the pt-slope of a' is the same as that of a" in T' . 

Because pT' is not in the interior of a' a" and a' and a" are different points, when moving 
a point from a' to a" on a'a" of the pt-slope of the moving point in T' changes 
monotonically (e.g., nondecreasing) . Recall that ^0(^2) B{T'), so ta(^2) B{dT'). 
So there must exist a point a* on a'a" such that the pt-slope of a* is larger than the 
minimum of those of a' and a" in T'. But, this fact and that the pt-slope of the moving 
point along a'a" changes monotonically together contradict with that the pt-slope of a' 
is equal to that of a" in T' . Thus, this possibility cannot occur. 

(b) Only one endpoint of ta{C2) is on dT' . The analysis here utilizes some analysis techniques 
for Case 1 (specifically for the subcase of a = g S yLZL\{zL})- As in Case 1, we will show 
that the pt-slope of l{p) in R' at the moment ^1 is smaller than the pt-slope of l{p) in T' 
at the moment ^2- Since l{p) is the directed tangent line of both the pseudo-triangles R' 
and T' at p, by Observation HI the pt-slope of l{p) in R' (at the moment ^1) must be no 
smaller than the pt-slope of l{p) in T' (at the moment ^2), which incurs a contradiction. 
Let /(ta(^2)) be the directed line that contains ta(^2) with the same direction as ta{C2)- 
We consider l{ta{^2)) as a physically directed line that is not associated with any time 
moment. Recall that ta(?2) and t'^ have opposite directions, l{ta{£,2)) and t'^ have opposite 
directions. We claim that the /(ta(^2))-slope of l{p) is larger than the ^R'-slope of l{p) at 
the moment ^1. Indeed, at the moment ^1, the 6/j'-slope of t'^ is the pt-slope of the point 
a on dR', which is larger than zero and less than it. Since t'^ and l{ta{C2)) have opposite 
directions, the ^(ta(^2))-slope of bji' is less than vr and larger than zero. Note that l{p) 
is the directed tangent line of R'. So the -slope of l{p) is the pt-slope of l{p) in R' (at 
the moment ^1), which is less than vr. Recall that ta{(,2) crosses l(p) from left to right, 
implying the /(ta(^2))-slope of l{p) is less than vr. Therefore, the Z(ta(^2))-slope of l{p) is 
the sum of the l{ta{S,2))-slope of bjii and the ft^j'-slope of l{p). Since the /(ta(^2))-slope 
of bn' is larger than zero, the claim is true. Since the bji'-slope of l{p) is the pt-slope of 
l{p) in R', we obtain that the /(ta(^2))-slope of l{p) is larger than the pt-slope of l{p) in 
R' at the moment ^1. 

Let a' be the endpoint of ta{(,2) on dT'. Since the direction of l{ta{(,2)) is the same as 
ta{£,2), l'i^aiC2)) is the directed tangent line of T' at a'. Recall that fa(^2) crosses l{p) 
from left to right, implying the ^(ta('?2))-slope of l{p) is less than vr. By Observation U 
the Z(ia(^2))-slope of l{p) is no bigger than the pt-slope of l{p) in T' (at the moment ^2)- 
We thus obtain that the pt-slope of l{p) in R' at the moment ^1 is smaller than the 
pt-slope of l{p) in T' at the moment ^2- Consequently, this possibility cannot occur. 

(c) Neither endpoint of ta{(,2) is on dT'. In this situation, since ta{(,2) intersects l{p) (at a) 
in the interior of T', ta{(,2) must cross dT' somewhere. Note that l{p) intersects ta(^2) (at 
a) before q and ta(C2) crosses l{p) from left to right. Recall in our proof by contradiction 
we assume the atom A is awake on dT' and the point p £ A lies on Awake[T'] = Sffior' 
(see Fig. [T5]) . Note that gj" always lies on xjvzt' and so does the point p. Since the 
interior of T' is free of obstacles and ta{^,2) crosses l{p) from left to right, similar to the 
analysis for Case 1, there must be a directed free bitangent t'^ G B{dT') lying on xjvp 
or pzT' such that ^0(^2) crosses t'^ from left to right (^0(^2) may also cross zjyxr', but 
we are not interested in that). Thus t'^ crosses ta(^2) from right to left. However, since 
t'l G B{T') and ta(^2) B{T') for the good pseudo-triangulation T' at the moment ^2, 
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Figure 15: Illustrating the situation when no endpoint of ta{^2) is on dT': ta{(,2) crosses l{p) from left to right. 

by the third property of good pseudo-triangulation, t[ should cross ta{£,2) from left to 
right. But this is a contradiction. Hence, this possibility cannot occur. 

We conclude that the case when q lies on a bitangent in B{dT') cannot occur. 

Below we discuss the case when q lies on an obstacle. For simplicity, here we view qas a. physical 
point not associated with any time moment. Consider the position of q with respect to the chain 
y^fzji' on dR'. There are two subcases to consider: (i) q lies on yRf^'] (ii) q does not lie on 
(i.e, l{p) crosses a free bitangent in yuf^')- We show below that neither subcase can occur. 

(i) q lies on yjfzji'. Let l'{q) be the directed tangent line of R' at q (at the moment .^i). Recall 
that l{q) is the directed tangent line of T' at q (at the moment ^2) and l{q) crosses l{p) from 
left to right. Since q lies on an obstacle, l{q) and l'{q) both lie on the same undirected line 
but may have opposite directions. The following analysis is very similar to that for Case 1 
(specifically for the subcase oi a = q G Vlzl \ {^l})- 

If the point q has not been reversed since the moment ^1, then l{q) has the same direction 
as l'{q). As analyzed before, the pt-slope of l{p) is smaller than that of l'{q) in i?', and thus 
l{p) crosses l\q) from left to right, or equivalently, l'{q) crosses l{p) from right to left, which 
contradicts with that l{q) (= l'{q)) crosses l{p) from left to right. 

If the point q is reversed (at most once by Lemma [6|) during the time from ^1 to ^2, then l'{q) 
and l{q) have opposite directions. By the same analysis as for subcase (iv (b)) of the former 
case (i.e., the case when q lies on a bitangent in B{dT')), we can show that the pt-slope of 
l{p) in R' at the moment ^1 is smaller than the pt-slope of l{p) in T' at the moment ^2) which 
contradicts with Observation [H We omit the details. Thus, this subcase cannot occur. 

(ii) q does not lie on yifz^i (i.e, l{p) crosses a free bitangent on yn/^'). Clearly, l{p) must cross 
a free bitangent on yj^zji' before arriving at q. Then, the analysis follows in exactly the same 
way as that for subcase (iv) of the former case, and we can conclude that this subcase cannot 
occur either. 

In summary, we prove that the atom A cannot be awake at the moment ^2 of the algorithm. 
For Case 2.2, we conclude that when T' = Rtri{bT'), the atom A on dT' cannot be dequeued 
due to the flip of bx' ■ 

Based on the above detailed case analysis. Sub-lemma [TOT a) is proved. 

By an analogous analysis, we can also prove the following statement, called Sub-lemmallOlfbl. 

Sub-lemma [TOj fb). Suppose A is an atom on dT of a pseudo-triangle T with T = Ltri(bT) in a 
good pseudo-triangulation T and A is dequeued due to a flip operation on bx- Also, suppose at any 
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later moment before the reversal of A, A lies on dT' of a pseudo-triangle T' with T' = Ltri{bT') in 
another good pseudo-triangulation T'. Then A cannot be dequeued due to the flip of bx'- 

By combining S ub- lemmas llOlf a) andliorb). we conclude that for any dequeued atom, it can be 
dequeued at most twice before its next reversal (if any). An explanation of this is as follows. 

Suppose at the moment ^i, an atom A on dTi of a pseudo-triangle Ti in a good pseudo- 
triangulation 71 is dequeued for the first time due to the flip of 6ti- We assume Ti = Rtri{bTi) (the 
case for Ti = Ltri{bTi) can be analyzed similarly). If A will not be dequeued again in the algorithm, 
then we are done. Otherwise, suppose at a later moment ^2 (A has not been reversed), A is on 
dT2 of a pseudo-triangle T2 in a good pseudo-triangulation Ti and A is dequeued for the second 
time due to the flip of 6t2- By Sub-lemma fTOT a). the case T2 = Rtri{bT2) cannot occur. Thus, only 
T2 = Ltri{bT2) is possible. For any moment ^3 after suppose A is on 5T3 of a pseudo-triangle 

in a good pseudo-triangulation Ts and A has not been reversed. Now, if = RtriibT^), then 
by Sub- lemma [TOT a) . A cannot be dequeued due to the flip of br^- But if T3 = Ltri{bT3)-, then by 
Sub- lemma [TOlf bl . A cannot be dequeued either. 

Thus, A can be dequeued at most twice in Phase I of the entire algorithm before its next reversal 
(if any). Lemma [TO] then follows. 

Note: It appears possible to show that there is at most one dequeue per atom before its reversal. 
But, proving this seems to make the already long and complicated proof of Lemma [TOl even longer 
and more complicated, and this stronger statement, although nicer, is not essential to our result. 

In addition, we briefly discuss why R' (= Rtri(b*)) is also Ltri{b'^) in Case 2 with the help of 
Lemma [TTl In Case 2.1 (i.e., xr = Head{b'j^)), this is obviously true. We discuss Case 2.2 (i.e., 
xji = Head{b)) below. It suffices to show that b'j^ lies on xj^yfi of dR. Assume to the contrary 
that b'j^ does not lie on XRyR of dR. Then, xj^ijr must be an obstacle arc on an obstacle, say P. 
Thus, it is easy to see that b'j^ must be tangent to P and the tangent point on P is i/r, and b'j^ lies 
on yRZR. Further, since in Case 2, p* does not lie on the obstacle arc between b and 6'^, we have 
P* XRyR, and thus b'^ G VrP* ■ Note that vrp* is part of xrOr- Therefore, b'^ S VrMr- However, 
since yR G xrcir, by Lemma [TTl the portion yRqR of dR is an obstacle arc, contradicting with 
^'r £ yR<lk since b'^ is a free bitangent. Hence, b'^ must lie on XRyR and R' is Ltri{b'^). 

B.3 Bounding the Number of Enqueue Operations in Phase II 

In this section, we prove ue < nq + nD + ns + k and ng = 0{n + k). Consequently, due to 
nq = 0{n + k), flu = 0{n + k), and k = we obtain ue = 0{n + k) and Lemma|4]thus follows. 

Let Q be the set of all reversed atoms in the entire algorithm, D be the set of all dequeue 
operations in Phase I, and S be the set of special enqueue operations in Phase II that will be 
defined later. Thus, nq = \Q\, no = \D\, and ns = \S\. 

To prove ue < nq + n£) + ns + k, we will show that every enqueue operation in Phase II 
corresponds to an element in D, Q, S, or B, and we charge the enqueue operation to that element; 
each such element will be charged only 0(1) times in the entire algorithm. Further, we will prove 
that ns = \S\ = 0{n + k). We discuss the three main cases individually. 

For Case 1, refer to the pseudocode Algorithm [1] in Section |Al There are three enqueue se- 
quences, i.e., LineslHllll andO and their corresponding charges are already shown in the pseudocode. 
We briefly explain why we can charge them in those ways. For LineO it is easy to see that ^lPr/ 
(recall pRi = Tail{b'^) is reversed due to the flip of b. So we can charge the enqueue on 'zlPr; to 
Q. For Line [H since b* S B, we can charge the enqueue on b* to B. For Line [HI recall that when 
computing b* in Phase I, all atoms in q'p^p* have been dequeued (from AwakeMax[R]). Note that 
the atoms in the enqueued portion of Line [8] are all in q'^p*, implying that they have been dequeued 
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in Phase I. So we can charge them to D. 

For Case 2, again, we discuss the two subcases Case 2.1 (i.e., xr = Head{b'j^), see Fig. [8]) and 
Case 2.2 (i.e., xr = Head{b), see Fig. (H). 

For Case 2.1, refer to the pseudocode Algorithmic There are three enqueue sequences, i.e.. 
Lines [3l HI and [H and their corresponding charges are already shown in the pseudocode. For Line 
[3l note that zlxr is reversed due to the flip of b, so we charge the enqueue operations on zlxr 
to Q. For LineHl we charge the enqueue on b* to B. For Line El note that the atoms in yup* are 
dequeued when computing b* in Phase I, so we charge the enqueue to D. 

For Case 2.2, refer to the pseudocode Algorithm [3 There are five enqueue sequences, i.e., 
Lines O [3 [TTl [151 and [19l For Line [3l which is for the case when qRi G Pr'P* , note that the 
enqueued portion is XRiqRi. Since qri G Pr'P*, Qr' is before p*, and thus all the enqueued atoms 
in Line [3] are dequeued for computing b* in Phase I. Therefore, we can charge them to D. Line 
[7] is trivial. For Line [TTl again, the atoms in the enqueued portion yRp* have been dequeued for 
computing b*, we charge those enqueue operations to D. For Line[15l note that pr' = Head{b'^. 
Since xr = Head{b), the enqueued portion xrPri is reversed due to the flip of 6, so we can charge 
the enqueue to Q. The enqueue in Line [191 needs special treatment. Below we discuss that the 
enqueued atoms there have some special properties and we call those enqueue the special enqueue 
operations and charge them to S. Later, we will prove IS*! = 0{n + k). 

We assume that the reader has read the detailed algorithm discussion for Algorithm [3] in Section 
[Al Note that Line HH is for the case wl = ql- First, note that the enqueue sequence in Line [191 
is on WRiqi of dR' , which is part of qR/qZ regardless of whether wri is yRi or qRi. We claim that 
WR^qi is part of yLQL on dL. To prove this claim, it suffices to show that wri is after yi on dL. 
Clearly, wri is always after yRi. Recall that yRi is q* or y^. If yR' is q*, then q* is after yL on dL; 
else, yR' is yL. In either case, yR' is always after yi on dL. Consequently, wr' is after yi on dL. 
The claim thus follows. Hence, the special enqueue sequence is on WR/qL, which is part of both 
qRf^ and y[qZ. 

Due to Wl = qL, yL is before qL on dL, and in other words yL lies on xlQl- By Lemma [TTl VlQl 
is an obstacle arc and thus WR'qL also lies on that obstacle. For any point p on an atom A of WR/qL 
on dR', let l{p) be the directed tangent line of the pseudo-triangle R' at p. Suppose we move from p 
along l{p) towards its inverse direction (resp., the direction of l{p)), and let a (resp., a') be the first 
point encountered on any obstacle in V; we call the point a (resp., a') the backward V-view (resp., 
forward V-view) of p. Let Vbacki^) (resp., VforiA)) denote the set of backward (resp., forward) 
"P-view points of the points of A. Since xrPr' (recall pRi = Head{b'j^)) on dR is reversed due to the 
flip of 6, by Observation [H xrPr' is an obstacle arc. An easy but critical observation is that for any 
point on qR/qZ, its backward "P-view is on xrPr' (see Fig. [9]). Since the enqueued portion WR/qL is 
part of 'qRiqZ-, for any point on WRiqL \ {qr'}, its backward P-view is on xrPr>, which is reversed 
due to the flip of b. To summarize what have been deduced above, we have (i) the portion of dR' 
involved in the special enqueue sequence is WR/qZ, which is an obstacle arc, and (ii) the backward 
P-view points of all points on WRiql \ {qr'} lie on xrPri, which is an obstacle arc reversed due to 
the flip of b. We then have the following observation. 

Observation 5 Consider a flip operation on a minimal bitangent b. If an atom A on dR' with 
R' = Rtri{(p{b)) is involved in a special enqueue operation (for processing the flip ofb), then A is 
an obstacle arc and Vhacki^) is an obstacle arc reversed due to the flip ofb. 

In addition, if an atom A is an obstacle arc, say, on the obstacle P, then for any point p G A, 
the tangent line of P at p is also the tangent line of the pseudo-triangle (at p) on which A lies, 
and vice versa. Thus, as long as the atom A is not reversed in the algorithm, both Vbacki^) and 
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VforiA) will not change. After A is reversed (if this ever happens), VbackiA) and Vfor{A) still refer 
to the same two obstacle arcs but switch names with each other. 

The bound of |5| will be discussed later in Lemma [T2l Some special enqueue operations also 
appear in Case 3. 

For Case 3 (see Fig. fTO]) . refer to the pseudocode Algorithm HI There are three enqueue se- 
quences, i.e.. Lines and [HI For LineO since xrp* is reversed due to the flip of 6, we charge the 
enqueue to Q. Line H] is trivial. For Line El we explain below that the enqueue can also be viewed 
as special enqueue and charged to S. 

Note that Line [8] is on the case wl = ql. Thus, i/l £ xiqi- By Lemma [HI yi/lL is an obstacle 
arc. Recall that in Case 3 xr' is either i/l or q*. If xji' = q* , q* is on yLzZ', otherwise, xri = ul- In 
either case, the enqueued portion in Line[8l i.e., xj^/qZ, lies on yiQl, which is an obstacle arc. Note 
that the backward V-view of any point on xj^/qZ is on xjip* (xji = Head{b), see Fig. [TOj) . which is 
reversed due to the flip of b and is an obstacle arc by Observation [H 

In summary, we have: (i) the enqueued portion xjiiqi in Line [8] is an obstacle arc, and (ii) the 
backward P-view points of all points on xj^/qZ lie on xrp*, which is an obstacle arc reversed due 
to the flip of b. Thus, Observation [5] also applies to the enqueue on xfj^^qZ in Line[8j Therefore, we 
also treat the enqueue as special enqueue and charge them to S. 

We have finished the discussion on charging the enqueue operations in Phase II for all cases. 
It remains to show l^l = 0{n + k). For this, we first give a similar observation on the (possible) 
special enqueue operations on dL' with L' = Ltri{ip{b)) for the flip operation on b. 

Observation 6 Consider a flip operation on a minimal bitangent b. If an atom A on dL' with 
L' = Ltri(ip{b)) is involved in a special enqueue operation (for processing the flip ofb), then A is 
an obstacle arc and VforiA) is an obstacle arc reversed due to the flip ofb. 

Based on Lemma [6] and Observations [5] and [H we prove the following lemma. 

Lemma 12 The number \S\ of all special enqueue operations in Phase II of the entire algorithm 
is 0{n + k). 

Proof: Consider an arbitrary atom ^ in a good pseudo-triangulation at a moment during the 
algorithm {A may have been reversed). First, it is important to note that A can be involved in a spe- 
cial enqueue operation only if A is an obstacle arc and the special enqueue operation is for processing 
a flip operation on a minimal bitangent b such that A is on dT with T G {Rtri{ip{b)) , Ltri{ip{b))} . 
In the following, we show that the atom A can be involved in at most two special enqueue opera- 
tions before its next reversal (if any). If A is not an obstacle arc, the above statement simply holds. 
We assume A is an obstacle arc. We assume that the discussion below is on the time period after 
the moment and before the next reversal of A (if any), unless otherwise stated. 

If the atom A is not involved in any special enqueue operation in the algorithm later on, then 
we are done. Otherwise, suppose at a later moment > ^q, A is involved in a special enqueue 
operation for processing the flip of a minimal bitangent ti. Without loss of generality, we assume 
that A is on dRtri{ip{ti)). Then, by Observation [5l the obstacle arc Vback{A) is reversed due to 
the flip of ti. By Lemma [6l Vhacki^) cannot be reversed again later in the algorithm. 

If A is not involved in any special enqueue operation in the algorithm after the moment ^i, then 
we are done. Otherwise, suppose at a later moment ^2 > Ci) ^ is involved in a special enqueue 
operation for processing the flip of a minimal bitangent t2- Because A has not been reversed since 
the moment as discussed earlier, both VbackiA) and Vfori^) do not change. Then, A cannot 
be on dRtri{ip{t2)) since otherwise, by Observation [H the obstacle arc Vback{A) would be reversed 
again due to the flip of t2- Hence, A can only be on dLtri{ip{t2)). By Observation [6l the obstacle 
arc VforiA) is reversed due to the flip of t2- By Lemma [Gj VforiA) cannot be reversed again later. 
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Consider any flip operation on a minimal bitangent ts in the algorithm at any later moment 
^3 > ^2- Because A has not been reversed since the moment (,o, both Vbacki^) and Vfor^A) do 
not change. No matter whether A is on dRtri{(p{t3)) or dLtri{ip{t3)), A cannot be involved in any 
special enqueue operation for processing the flip of ts, since otherwise, by Observations [5] and El 
the obstacle arc VtackiA) or VforiA) would be reversed again due to the flip of ^3. 

Therefore, we obtain that A can be involved in at most two special enqueue operations before its 
next reversal (if any). By Lemma [6l A can be reversed at most once. We now claim that A can be 
involved in at most two special enqueue operations in the entire algorithm. Indeed, if A is involved 
in two special enqueue operations before its reversal (if any), then both VbackiA) and Vfor{A) are 
reversed before the reversal of A. Note that after A is reversed, Vbacki^) and VforiA) refer to the 
same two obstacle arcs but switch names with each other. Since neither of Vtacki^) and Vfori-A) 
can be reversed again, the reversed atom A cannot be involved in any special enqueue operation in 
the rest of the algorithm. If A is involved in one special enqueue operation before its reversal, then 
exactly one of Vbacki^) and Vfor{A) is reversed before the reversal of A. After A is reversed, only 
one of Vback{A) and Vfori^) (with their names switched with each other) can possibly be reversed, 
and thus in this situation the reversed A can be involved in at most one special enqueue operation 
in the rest of the algorithm. Finally, if A has not been involved in any special enqueue operation 
before its reversal, then the claim simply holds. Therefore, the claim is true and the atom A can 
be involved in at most two special enqueue operations in the entire algorithm. 

Since the number of all atoms in the algorithm is 0(n + k), the number of all special enqueue 
operations in Phase II of the entire algorithm is 0{n + k), i.e., |5| = 0{n + k). □ 

C The Bounded Degree Decomposition of the Free Space 

Recall that S = {Si, . . . , Sh} is a set of h splinegons of totally n vertices with pairwise disjoint 
interior. In this section, we compute a bounded degree decomposition BDD[F) of the free space 
In the following, we first define BDD{T) and then present our algorithm for it. 

C.l Defining a Bounded Degree Decomposition 

As preprocessing, we perform a monotone cut on the edges of the splinegons in 5, as follows. For 
each splinegon edge e, if one or both of its topmost and bottommost points lie in the interior of 
e, then we add these points (at most two) as new splinegon vertices, which divide the original 
edge e into several new edges (at most three). Since each splinegon edge is of 0(1) complexity, 
this monotone cut can be done in 0{n) time. After the cut, S contains at most 3n vertices. For 
convenience, with a little abuse of notation, we still use n to denote the number of vertices of S after 
the monotone cut. From now on, we assume that the monotone cut has been done on all splinegons 
in S. Recall that the splinegons in S have pairwise disjoint interior; for ease of exposition, we 
assume that they also have pairwise disjoint boundaries. For any object A in the plane, let dA 
denote its boundary. 

Recall that 7^ is a large rectangle containing all splinegons in S. Let dJ- denote the boundary 
of J-, i.e., the union of the boundaries of the splinegons in S and IZ. A diagonal is an open line 
segment in the interior of J- with its two endpoints on dJ-. A hounded degree decomposition of the 
free space J-", denoted by BDD(J^), is a decomposition of J- into 0(n) bounded degree regions (or 
simply regions) each with at most four sides and with at most three neighboring regions by adding 
0{n) non-intersecting diagonals (see Fig. [3]). Two regions are neighboring if they share a diagonal 
on their boundaries. Each region has at most four sides and each side is either a diagonal or (part 
of) a splinegon edge. Thus the complexity of each region is 0(1) (that is why we call it a "bounded 
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Figure 16: Illustrating the horizontal visibility map of Figure 17: Illustrating the decomposition of a trapezoid: 
a splinegon. (a) The black points on the base cd are vertices; (b) a 

decomposition of the trapezoid. 



degree region"). BDD{T) serves the same purpose as a triangulation in the polygonal domain case. 

For a single simple splinegon S, a linear time algorithm was given in [27] for computing a 
bounded degree decomposition of S. 

C.2 Computing a Bounded Degree Decomposition 

Our algorithm for computing BDD{T) consists of two main steps. First, we find h non-crossing 
diagonals to connect all splinegons in S and IZ together to form a single simple splinegon S* . 
Second, we apply the algorithm in [27] to computing a bounded degree decomposition of 5*, which 
is BDD{T). Below, for completeness and easy understanding of our algorithm, we briefly discuss 
the linear time decomposition algorithm in [27] for a single splinegon. 

Let 5 be a splinegon of n vertices. The algorithm in [27j first decomposes S into a set of 
horizontal trapezoids by computing the horizontal visibility map and then further decomposes each 
trapezoid into (bounded degree) regions as the final decomposition of S. Below are some details 
of it. Again, the topmost and bottommost points of each splinegon edge are treated as vertices of 
S. Clearly, there are 0{n) vertices on dS. As for the simple polygon case, the horizontal visibility 
map is to draw a horizontal line segment through each vertex of 5, extending the segment so long 
as it does not properly cross dS (see Fig. [TBI) . The visibility map of S adds 0{n) new vertices on 
dS and divides S into 0{n) trapezoids (each with curved sides and line segment bases). As shown 
in [27] . Chazelle's algorithm [3] can be used to compute the visibility map on the splinegon S in 
0{n) time. Since there may be multiple vertices in the interior of a trapezoid base, a trapezoid may 
have many neighbors. The next step is to further decompose each trapezoid into (bounded degree) 
regions such that each region has at most three neighbors. There are many ways to decompose a 
trapezoid into (bounded degree) regions. In [27], one algorithmic approach for this task was given 
and some of the cases were discussed. For example, consider a trapezoid abed with ab and cd as 
bases and ad and be as (curved) sides, as in Fig. [TW a). Suppose there are multiple vertices in the 
interior of the base ed. A further decomposition of the trapezoid is shown in Fig. [TTt^ b). Refer to 
[27] for more details. As can be seen from the example, the algorithm in |27j tends to make more 
regions with three sides. We do not have to do this although doing so as in [27] does not affect the 
performance of our algorithm asymptotically. Actually, a region may be a four-side trapezoid, but 
it has at most three neighbors. The key is to ensure that each region has at most three neighbors 
and is of 0(1) complexity. 

We then discuss the first step of our decomposition algorithm on J^, i.e., finding h non-crossing 
diagonals to connect all splinegons in S and TZ together into a single simple splinegon. Our approach 
generalizes the triangulation algorithm in [2] for the polygonal domain case. We first define a 
visibility tree for S, as follows. For each splinegon Si G S, pick a point on its boundary (not 
necessarily a vertex) and draw a ray to the right until it hits either another splinegon in S or TZ. 
As shown in [2j, if we choose the origins of the rays carefully and view each splinegon as a node, 
we can then ensure that the resulting planar graph is connected and acyclic; actually, the resulting 
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planar graph is a visibility tree for S, denoted by Tms{<S). Clearly, Tms{<S) connects all splinegons 
of S and TZ into a single simple splinegon. Our task below is to compute Tyis{S). Note that in the 
special case when all splinegons in S are convex, a visibility tree Tyis{S) can be easily computed in 
0{n + hlogh) time by a sweeping algorithm. The general case is handled by the next lemma. 

Lemma 13 A visibility tree Tms{S) of S andlZ can he computed in 0{n + hlog^'^'^ h) time. 

Proof: Our approach generalizes the algorithm in [5], which computes a visibility tree in 0(n + 
/ilog^^*^ h) time for a set of h pairwise disjoint polygons of totally n vertices. 

To generalize the algorithm in [2j, we need to make sure that each of its components can be 
generalized. First, the algorithm in |2j makes use of the linear time algorithm in [17 \ for sorting 
the intersections (by their x-coordinates) of a horizontal line and an oriented Jordan curve. For 
any splinegon Si G 5, consider the problem of sorting the intersections of dSi with a horizontal 
line /. Due to the monotone cut, each edge of Si has at most one intersection with / which can be 
computed in 0(1) time. Further, the line / breaks up the boundary of Si into disjoint arcs either 
entirely below or above / such that the arcs above (resp., below) the line still form a parenthesis 
system, as shown in [2l[T7|. Hence, the algorithmic scheme in [IT] and the time analysis are still 
applicable to our problem. In summary, the intersections of dSi and the line / can be computed in 
linear time (in terms of the number of vertices of Si); let m be the number of such intersections. 
Then these intersection points on / can be sorted in 0{m) time. 

With the above sorting algorithm in hand, the following more general sorting problem can be 
solved by using the approach in [2]: Given a subset of h' {h! < h) splinegons in S with totally n' 
(n' < n) vertices and a horizontal line /, suppose the total number of vertices in these splinegons 
is n'; the goal is to sort the intersections of I with these h' splinegons. All intersections can be 
computed in 0{n') time. Let m be the number of such intersections. Then by following the 
algorithmic scheme in and using our sorting procedure for a single splinegon case, these m 
intersections can be sorted in 0{m + h' log hi) time. 

In addition, the algorithm in [2] needs a point location data structure |101 [23] on a simple 
polygon, constructed in linear time, for answering each point location query in logarithmic time. In 
our problem, correspondingly, we need such a point location data structure on a simple splinegon. 
As shown in [101 [27] , the data structure in [lOj can be made work on a simple splinegon with the 
same performance as for the simple polygon case. 

It is also easy to verify that other parts of the algorithm in [2] are all applicable to our problem. 
Therefore, the visibility tree Tms{S) for S can be computed in 0{n + /ilog^"*"^ h) time. □ 

By Lemma [13] and the linear time decomposition algorithm for a simple splinegon [27] , the 
following result follows. 

Theorem 2 A hounded degree decomposition of the free space T among the splinegons of S can he 
computed in 0{n + hlog^'^'' h) time. If all splinegons in S are convex, then the decomposition can 
be computed in 0{n + h log h) time. 

D Computing the Corridor Structure 

In this section, based on BDD{T), i.e, the bounded degree decomposition of the free space T, we 
compute a corridor structure that helps reduce our original problem on SPSD to an instance of the 
convex SPSD. Our algorithm generalizes the approach in [21| for the polygonal domain case. 

Recall that both s and t are considered as two special splinegons in S. In addition to the 
splinegon vertices, the endpoints of the diagonals of BDD{T) are also treated as the vertices of 
BDD{F). Note that BDD{F) has 0{n) vertices. Let G{F) denote the planar dual graph of 
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BDD{F), i.e., each node of G{F) corresponds to a region in BDD{F) and each edge connects two 
nodes of G{J-) corresponding to two regions sharing a diagonal. Because the sphnegons in S are 
pairwise disjoint, the dual graph G{J-) is clearly connected, and an s-t path among the splinegons of 
S always exists. Since BDD[T) is a planar structure and each region in BDD{T) has at most three 
neighbors, G{J-) is a planar graph whose vertex degrees are at most three. As in the polygonal 
domain case, at least one node dual to a region incident to each of s and t is of degree three. 

Based on G{F), we compute a planar 3-regular graph, denoted by G^ (the degree of each node 
in it is three), possibly with loops and multi-edges, as follows. First, we remove every degree-one 
node from G{J-) along with its incident edge; repeat this process until no degree-one node exists. 
Second, remove every degree-two node from G[F) and replace its two incident edges by a single 
edge. The resulting graph is G^ (e.g., see Fig.[3j). By a similar argument as in [21] for the polygonal 
domain case, we can show that the resulting G^ has h + 1 faces, 2h — 2 nodes, and 3/i — 3 arcs. Each 
node of G^ corresponds to a region of BDD{F), which is called a junction region (e.g., see Fig. [3]). 
Removal of all junction regions from G'^ results in 0{h) corridors^ each of which corresponds to one 
edge of G^. 

The boundary of a corridor G consists of four parts (see Fig. S]): (1) A boundary portion of a 
splinegon Si G 5, from a point a to a point 6; (2) a diagonal of a junction region from 6 to a point c 
of a splinegon 5*2 G 5 (it is possible that Si = S2); (3) a boundary portion of the splinegon 5*2 from 
c to a point d] (4) a diagonal of a junction region from d to a. The two diagonals be and ad are called 
the doors of C. Note that the corridor C itself is a simple splinegon. Let \C\ denote the number of 
vertices of BDD{T) on dC. Note that a shortest path between two points inside a simple splinegon 
can be found in linear time [27j. Therefore, in 0(|C|) time, we can compute the shortest path 7r(a, h) 
(resp., 7r(c, d)) from a to 6 (resp., c to d) inside G. The region He bounded by 7r(a, 6), 7r(c, d), and 
the two diagonals he and da is called an hourglass, which is open if 7r(a, h) n 7r(c, c?) = and closed 
otherwise (see Fig. [3|). If He is open, then both 7r(a, h) and 7r(c, d) are convex and they are called 
the sides of He] otherwise, He consists of two "funnels" and a path vrc = vr(o, b) R 7r(c, d) joining 
the two apices of the two funnels, called the corridor path of G. Each funnel side is also convex. 
We process all corridors as above. The running time for processing all corridors is linear in terms 
of the total number of vertices of all corridors, which is at most the number of vertices of BDD{T), 
i.e., 0{n). Therefore, the running time for processing all corridors is 0{n). 

Let Q be the union of all junction regions and hourglasses. Then Q consists of 0{h) junction 
regions, open hourglasses, funnels, and corridor paths. Let 7r(s, t) be a shortest s-t path for the 
original problem SPSD. As shown in [2lJ, Ti{s,t) must be contained in Q. Consider a corridor G. 
If 7r(s, t) contains an interior point of G and neither s nor t lies on dG, then the path 7r(s, t) must 
cross both doors of C, i.e., it enters G from one door and leaves G from the other. Further, if the 
hourglass He for G is closed, then the corridor path of G must be contained in 7r(s,t). When He 
is open, since both sides of He are convex with respect to the interior of He-, if vr(s,i) intersects 
both sides of He, then it must contain a common tangent of the two sides such that vr(s,t) goes 
from one side of He to the other side via that common tangent. 

With all the properties above, let Q' be the set of the sides of all open hourglasses and the sides 
of all funnels in Q. Then Q' consists of 0{h) convex chains. As in the polygonal domain case |21j . 
by somehow extending each convex chain in Q' , we can partition the space into a set S' of 0{h) 
convex splinegons of totally 0(n) vertices with pairwise disjoint interior, such that a shortest path 
7r(s,t) for the original SPSD is also a shortest s-t path that does not intersect the interior of the 
splinegons in S' but possibly contains some corridor paths. Therefore, other than the 0{h) corridor 
paths, we have reduced our original SPSD problem to an instance of the convex SPSD. As in |21j . 
the convex splinegon set S' can be easily obtained, say in 0{n + hlogh) time. 
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E Computing a Shortest Path in the Relevant Visibihty Graph 



In this section, we compute a shortest path from s to i in G, which is the relevant visibihty graph 
of the 0{h) pairwise disjoint convex sphnegons in S' with totally 0(n) vertices. Recall that k is 
the number of the free common tangents of all splinegons in S' . For convenience, we assume that 
the number of convex splinegons in S' is h and the total number of splinegon vertices is n. Let 
S' = {Si,S2, . . . , Sh}- 

To find a shortest path from s to t in the graph G, since G has 0{k) nodes and 0{k) edges, 
simply running Dijkstra's algorithm on G would take 0(A;logA;) time. To avoid the log A; factor, we 
transform G to a coalesced graph G'^ such that: (1) has only 0{h) nodes and 0{k) edges; (2) a 
shortest s-t path in G corresponds to a shortest s-t path in G^, which can be found in 0{h\ogh + k) 
time. This approach is quite similar to that in [5j for computing a shortest s-t path among n convex 
pseudodisks of 0(1) complexity each. In general, the approach in [5j relies only on the convexity 
of the objects involved and thus is applicable to our problem setting. Note that the idea of using 
a coalesced graph was first proposed by Hershberger and Guibas [15], but the definition of the 
coalesced graph and its construction in [5] are both different from those in [15j. We extend the 
method in [5] to solving our problem in the splinegon setting. 

As the approach in [5j, a key to our algorithm is to compute a set of 0{h) "distinguished points" 
on the boundaries of the splinegons in 5', which are then used to construct G'^. By a proof similar 
to that in [5J, a set of 0{h) distinguished points can be obtained easily once the Voronoi diagram 
of the convex splinegons in S' is available. Denote by VD{S') the Voronoi diagram of the h convex 
splinegons in S' . The next lemma follows from the results in [5] (see [5j for more details). 

Lemma 14 /5/ After the Voronoi diagram VD[S') is built, the coalesced graph G^ with 0{h) nodes 
and 0{k) edges can be constructed in 0{n + k + hlogh) time. 

It remains to describe how to compute VD{S'). In Appendix lE.H we will show that VD(S') 
can be computed in 0(n + hlogh) time, which is optimal. Thus, we have the following result. 

Theorem 3 A shortest s-t path for the convex SPSD can be found in 0{n + hlogh + A;) time, 
where k = 0{h'^) is the number of free common tangents among the convex splinegons of S' . 

It seems that there is another way |14j to compute the distinguished points although the time 
bound would remain unchanged. 

E.l The Voronoi Diagram of Convex Splinegons 

In this section, we compute the Voronoi diagram VD{S') for a set S' of h pairwise disjoint convex 
splinegons of totally n vertices. To our best knowledge, no efficient algorithm was given previously 
for it. By extending Fortune's sweeping algorithm [IT], one may obtain an 0(n + /i log /i log n) time 
solution for it. For the polygon case (i.e., all splinegons in S' are convex polygons), VD{S') can be 
computed in 0{n + hlogh) time |26j . We show that by generalizing the algorithm in [26] , VD{S') 
for the convex splinegon case can also be computed in 0{n + hlogh) time (as in j26j . we assume 
that the edges of each input splinegon are represented as a cyclically ordered list). 

In fact, as in [26j, we achieve a stronger result: The compact diagram (to be defined below) of 
the convex splinegons in S' can be computed in 0{h log n) time, from which VD{S') can be derived 
in an additional 0{n) time. Note that hlogn = 0{n + hlogh). As in [26] and to be discussed later, 
the compact diagram has several advantages over the "normal" Voronoi diagram. 

We first formally define the compact diagram of 5', denoted by GD{S'). We follow the termi- 
nology in [26j. Consider a convex splinegon S € 5', which is contained in a Voronoi cell of V7?(5'), 
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Figure 18: Illustrating the compact diagram of a set of convex splinegons: The black polygons are the cores of the 
splinegons. The Voronoi diagram of the splinegons is shown with dashed curves. 

say Cs- For each Voronoi vertex v on the boundary of Cs-, we draw a hne segment from v to 
its closest point, say py, on S. The segment vp^ is called the spoke from u to 5 and the point 
Ptj is called the spoke attachment point. If the cell Cs is unbounded, then there is a point on dS 
whose normal does not intersect the cell boundary; we view this normal as a spoke from an infinite 
Voronoi vertex to S. The core of 5 is the convex hull of all spoke attachment points on S. The 
compact diagram CD{S') is the union of all spokes and cores of the splinegons in S' (see Fig. [TH]) . 
Theorem m gives the algorithm for computing CD{S'). 

Besides its efficient construction, the compact diagram has several advantages over the Voronoi 
diagram. Note that while the Voronoi diagram VD{S') has 0{n) high degree (but still constant) 
algebraic curves, the compact diagram CD{S') consists of only 0{h) line segments. This feature 
makes the compact diagram easier and more efficient to display and represent. In addition, for 
applications in which knowing only two candidates for the closest splinegons is sufficient, the original 
splinegons can be discarded and only the 0{h) segments of CD{S') need to be stored, using 0{h) 
instead of 0(n) space. In |26) . two applications of the compact diagram were discussed, i.e., the 
post-office problem and the retraction motion planning problem, in which the sites were modeled 
as convex polygons. With our results, if the sites are modeled as convex splinegons, then the 
corresponding post-office problem and the retraction motion planning problem can be handled 
similarly with the same performance as in [26]. Our results may also find other applications. 

Theorem 4 The compact diagram CD{S') of the convex splinegons in S' can he computed in 
0{h\ogn) time, from which the Voronoi diagram VD{S') can he derived in an additional 0{n) 
time. 

Proof: We first focus on computing CD{S'). For this, we generalize the corresponding algorithm 
in |26j . Given a set V of h pairwise disjoint convex polygons of totally n vertices (each polygon is 
represented in a standard fashion), McAllister, Kirkpatrick, and Snoeyink gave an algorithm ^26] 
for computing the compact diagram of the convex polygons of V in 0(/ilogn) time. We refer to 
their algorithm as the MKS algorithm. We first sketch the MKS algorithm and then discuss our 
generalization of it on the convex splinegon set S'. 

Like Fortune's approach [11], the MKS algorithm is a sweeping algorithm, sweeping the convex 
polygons of V from (say) left to right. As in Fortune's algorithm, the Voronoi cell boundary 
maintained by the sweepline is called the sweep front or heach line, which consists of Voronoi 
edges between the sweepline and some polygons. A maximal connected portion of the sweep front 
between the sweepline and a single polygon is called a front arc. There are two types of events 
in the sweeping process. A site event occurs when the sweepline reaches the leftmost point of a 
polygon. A circle event occurs when the sweepline reaches the rightmost point of a circle that is 
tangent to three polygons of some consecutive front arcs on the sweep front. Clearly, there are 0{h) 
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site events and 0{h) circle events. The MKS algorithm focuses on computing the vertices of the 
Voronoi diagram as well as identifying the polygons that generate those Voronoi vertices. Two data 
structures are maintained by the sweeping algorithm: A balanced binary search tree that stores 
the sweep front and a priority queue that schedules the events in the order that the sweepline will 
encounter them. 

Generalizing Fortune's algorithm in a straightforward manner would take 0(/i.log/ilogn) time 
since every site event is processed in 0(log/ilogn) time. Specifically, at each site event, the 
sweepline is at the leftmost point, say p, of a polygon, and it needs to determine the front arc on 
the sweep front that is closest to p. A straightforward processing of this task takes 0(log/ilogn) 
time. Based on a critical observation [26J, the MKS algorithm handles each site event in O(logn) 
time. This observation states that the sweepline can be partitioned into disjoint intervals such that 
finding the nearest front arc to the point p is equivalent to locating in which interval the point p lies, 
which can be carried out in O(logn) time. Two subroutines heavily used in the MKS algorithm are 
spoke{p, A) and vertex{A, B, C). Given a convex polygon A and a point p outside A, spoke{p, A) 
returns the closest point on A to p, which can be implemented in O(logn) time by binary search. 
The subroutine vertex{A, B, C) takes three convex polygons A, B, and C as input and computes a 
finite or an infinite Voronoi vertex v such that the Voronoi cells for the polygons A, B, and C occur 
in a counterclockwise order around v. An advanced technique developed in [25], called tentative 
prune- and- search, is used to implement vertex{A, B,C) in O(logn) time. 

A generalization of the MKS algorithm to computing CD{S') for our convex splinegon set S' 
turns out to be quite natural. First, for each convex splinegon S of 5', we consider its leftmost 
and rightmost points as two new vertices of S, which may partition at most two splinegon edges 
of S into two new edges each. This step can be done in 0{h\ogn) time by binary search on each 
splinegon in S' . After this step, for each splinegon edge of the splinegons in S\ any vertical line 
can intersect it at most once. Next, we sweep the splinegons of S' from left to right. We define 
the site events and circle events similarly as in the MKS algorithm. Clearly, there are still 0{h) 
site events and 0{h) circle events. We also maintain the two data structures for storing the sweep 
front and for scheduling the events. Essentially, the MKS algorithm replies on the convexity of 
the objects involved. Since the splinegons in S' are convex, the scheme of the MKS algorithm is 
still applicable. For example, the critical observation used by the MKS algorithm for handling site 
events is based on the convexity of the polygons. In our problem, at each site event, the sweepline 
is at the leftmost point, say p, of a convex splinegon. To determine the nearest front arc on the 
sweep front that is closest to p, by following the same approach as for the MKS algorithm, we can 
also partition the sweepline into disjoint intervals such that finding the nearest front arc to the 
point p is equivalent to locating in which interval the point p lies. 

For the implementation details and the running time of our compact diagram algorithm, in 
general, since each splinegon edge is of 0(1) complexity, the operations on splinegons edges needed 
in the algorithm can be performed in the same order of time asymptotically as those on poly- 
gon edges in the MKS algorithm. Obviously, the subroutine spoke{p, A) can be implemented in 
O(logn) time by a binary search. For the subroutine vertex{A,B,C), the tentative prune-and- 
search technique [25] is also applicable to our problem. Specifically, this technique defines three 
continuous, monotone-decreasing functions that rely only on the convexity of the objects involved. 
In our problem, since all splinegons are convex, we can define three such functions in exactly the 
same way as those for convex polygons in the MKS algorithm. One basic operation needed in our 
algorithm is: For any point p on the boundary of a convex splinegon Si € 5', compute the normal 
of Si at p. Since each splinegon edge is of 0(1) complexity, this operation takes 0(1) time, as 
in the convex polygon case. Other operations can also be performed in the same order of time 
as their counterparts in the convex polygon case. We omit the details. Therefore, the subroutine 
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vertex{A,B,C) can be implemented in O(logn) time. 

We conclude that the compact diagram CD{S') can be computed in 0{h\ogn) time. 

Since each splinegon edge is of 0(1) complexity, as in [26], the Voronoi diagram VD{S') can be 
derived from CD{S') in an additional 0{n) time. The theorem thus follows. □ 

F Wrapping Things Up 

We now show how to find a shortest s-t path for our original SPSD problem on the splinegons set 
S. In Appendix [ni we build a corridor structure to obtain 0{h) corridor paths and a set S' of 0{h) 
convex splinegons of totally 0{n) vertices with pairwise disjoint interior such that a shortest s-t 
path for SPSD is also a shortest s-t path avoiding the convex splinegons in S' and possibly utilizing 
some corridor paths. In Appendix [Ej based on S' , we construct a coalesced graph such that: 
(1) has only 0{h) nodes and 0{k) edges; (2) a shortest s-t path avoiding the convex splinegons 
in S' corresponds to a shortest s-t path in C^. 

To compute a shortest s-t path for our original SPSD problem, our final step is to incorporate 
the 0{h) corridor paths into the graph to obtain a new graph such that a shortest s-t path 
for SPSD corresponds to a shortest s-t path in G^, as follows. 

Recall that a corridor path connects the two apices of two funnels. When building G*^, in 
addition to other distinguished points, we also treat all the 0{h) funnel apices as distinguished 
points. In this way, every funnel apex defines two vertices in (since every distinguished point 
defines two vertices in G'^). Consider a corridor path connecting two funnel apices u and v. Suppose 
the two vertices in G*^ defined by u (resp., v) are ui and U2 (resp., vi and ^2). After obtaining G'^, 
we add to G"^ eight directed edges e(iii, fi), e(iii, ^2), e(u2, wi), e{u2,V2), e{vi,ui), e{vi,U2), e(f2, ui), 
and e(f2, M2)) whose weights are the length of the corresponding corridor path. We do this for each 
corridor path, and then obtain GJj. Since there are 0{h) corridor paths, the graph GJ^, which still 
has 0{h) nodes and 0{k) edges, can be constructed in 0(n + /ilog h + k) time. A shortest s-t path 
for SPSD can then be computed by running Dijkstra's algorithm on G^, in 0{h\ogh + k) time. 

In summary, we have the following result. 

Theorem 5 Given a set S of h general splinegons of totally n vertices with pairwise disjoint 
interior in the plane, a shortest s-t path in the free space can he computed in 0{n + /ilog^^*^ h + k) 
time, where k is the size of the relevant visibility graph and is bounded by 0{h'^). 

G Proof of Lemma [1] 

Lemma[T] The initial good pseudo-triangulation To ofV can be constructed in 0{n + hlogh) time. 

Proof: To construct To of V, we modify and generalize the corresponding (0(n log n) time) al- 
gorithm in [31j for the set O. The algorithm is based on a rotational sweeping procedure, during 
which a visibility map, denoted by M{u') associated with the current rotational direction u' G [0, vr], 
is (implicitly) maintained. Consider a direction u'. Each obstacle Pi of V contains two extreme 
points each having a tangent line with slope u' such that Pi is between these two tangent lines. We 
denote by V{u') the set of such extreme points in all obstacles of V. 

We first define M(0) for u' = as follows (see Fig. [T9l) . For each extreme point in V{u'), we 
shoot one ray in the direction of u' and shoot another ray in the opposite direction until hitting 
some obstacles. The subdivision of the plane defined by all these rays and the obstacles of V is 
M(0), which can be viewed as similar to a trapezoidal decomposition of the free space. In general, 
for any u' > 0, M{u') is defined by the rays shooting from the points of V{u') in the direction of u' 
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Figure 19: The visibility map M(0). 



Figure 20: The visibiUty map M{n/2). 



and in the opposite direction until hitting some obstacles or bitangents of To obtained up to that 
moment of the rotational sweep, together with the obstacles of V. 

The algorithm first constructs M(0). For this, we compute the extreme point set V{0), which 
takes 0{n) time. Since the obstacles of V are pairwise disjoint, with a standard sweeping algorithm 
(from top to bottom), M(0) can be easily constructed in 0{n + hlogh) time. 

Starting at M(0), we rotate u' from to tt. During the rotation, the topology of M{u') is 
maintained implicitly. Specifically, the topology of M{u') does not change until u' becomes equal 
to the slope of a free bitangent & of To . When a new free bitangent b is detected, a "quadrangular" 
region (which contains two points of V{u')) in M{u') will disappear, and some rays shooting from 
V{u') will first hit h instead of some obstacles or other free bitangents of To already found (we 
then view these rays as hitting b without going further). At the same time, a "triangular" region 
(which contains only one point of V{u')) will emerge. If two triangular regions contain a same 
point of V{u') at their boundaries, then they are incident along a ray shooting from this point, and 
we merge these two regions by removing this ray. The resulting new visibility map is M{u') with 
the newly detected free bitangent 6 of To- We keep rotating u' in this manner. Fig. [20] illustrates 
M{u') for u' = 7r/2. Note that every new free bitangent is detected from two obstacles along the 
boundary of a quadrangular region. In this way, after the rotation of u' is over. To is obtained. 

The key to this procedure is to determine the rotation events, i.e., which free bitangent will 
be encountered next in the rotation. For this, the same strategy of the original (O(nlogn) time) 
algorithm for O in [31] is applied. The only difference is that a bitangent of two 0(1) complexity 
obstacles in [31] is found in 0(1) time, while in our problem we compute each free bitangent of To 
in O(logn) time. Since To has 0{h) free bitangents, computing all of them during the rotational 
sweep takes 0{h\ogn) time. Note that hlogn = 0{n + hlogh). 

In summary, constructing To takes 0{n + hlogh) time. The lemma thus follows. □ 



41 



